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The reaction systems, ^^Zn + ^**Ni, ^''Zn + ^^Mo, ^"Zn + ^^^Au, at 26A, 35A and 47A MeV, have 
been studied both in experiments with a A-k detector array, NIMROD, and with Antisymmetrized 
Molecular Dynamics model calculations employing effective interactions corresponding to soft and 
stiff equations of state (EOS). Direct experimental observables, such as multiplicity distributions, 
charge distributions, energy spectra and velocity spectra, have been compared in detail with those 
of the calculations and a reasonable agreement is obtained. The velocity distributions of a particles 
and fragments with Z > 3 show distinct differences in calculations with the soft EOS and the stiff 
EOS. The velocity distributions of a particle and Intermediate Mass Fragments (IMF's) are best 
described by the stiff EOS. Neither of the above direct observables nor the strength of the elliptic flow 
are sensitive to changes in the in-medium nucleon-nucleon (NN) cross sections. A detailed analysis 
of the central collision events calculated with the stiff EOS revealed that multifragmentation with 
cold fragment emission is a common feature predicted for all reactions studied here. A possible 
multifragmentation scenario is presented; after the preequilibrium emission ceases in the composite 
system, cold light fragments are formed in a hotter gas of nucleons and stay cold until the composite 
system underdoes multifragmentation. For reaction with ^^^Au at 47A MeV a significant radial 
expansion takes place. For reactions with ^*Ni and ®^Mo at 47A MeV semi-transparency becomes 
prominent. The differing reaction dynamics drastically change the kinematic characteristics of 
emitted fragments. This scenario gives consistent explanations for many existing experimental 
results in the Fermi energy domain. 

PACS numbers: 25.70Pq 

Keywords: Intermediate Heavy ion reactions, 47r detector array, antisymmetrized molecular dynamics model 
calculations, nuclear equation of state, in-medium nucleon-nucleon cross section, cold fragment emission, 
semi-transparency, radial expansion, elliptic flow 



I. INTRODUCTION 

One of the main aims of the study of heavy ion reac- 
tions is to explore the properties of nuclear matter at var- 
ious densities and temperatures. In intermediate heavy 
ion reactions (a few tens of MeV/nucleon to a few hun- 
dreds of MeV/nucleon), it is generally expected that the 
composite system of projectile and target nuclei is com- 
pressed and excited in the early stage of the reaction, and 
that the hot-dense nuclear system expands and breaks up 
by a multifragmentation process. Recently many studies 
have been undertaken to elucidate possible critical behav- 
iors for such matter U . In many of these studies, thermal 



and/or chemical equilibrium is assumed However in 
order to reach high enough in excitation energy and/or 
in density, the collisions become very violent and the col- 
lision processes become very complicated. Therefore it 
is indispensable to establish reliable microscopic dynam- 
ical models for the study of the properties of the highly 
excited matter produced in heavy ion reactions. 

In order to elucidate the reaction mechanism many mi- 
croscopic dynamical models for nuclear collisions have 
been proposed 0, Among such models, the molec- 
ular dynamics models are well suited to deal with the 
multifragmentation process. In the Classical Molecular 
Dynamics model, particles are treated as point particles 
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and their transport is governed by a classical equation of 
motion in a given mean field 0, 0, 01 ■ Nucleon-Nucleon 
collisions are taken into account and treated as hard- 
sphere scatterings without Pauli-blocking. In the Quan- 
tum Molecular Dynamics (QMD) model, each particle is 
described by a Gaussian wave packet 0, • Initial nu- 
clei are constructed by ensuring that there is less than 
one nucleon in each phase space cell of 1/h'^. During the 
time evolution of the wave packets the Pauli principle is 
respected only by the Liouvillc theorem of classical me- 
chanics. In the model nucleon-nucleon (NN) collisions are 
allowed and the Pauli blocking is treated in an approxi- 
mate manner. During the propagation of the wave pack- 
ets, however, the time evolution based on the classical 
equation of motion eventually leads the initial state into 
a Pauli-forbidden zone and the occupation number of nu- 
cleons in phase space often significantly exceeds 1/h^ 0- 
There have been several attempts to respect the Pauli 
principle more strictly during the propagation of the wave 
packets within the classical mechanics. Introduction of 
a Pauli potential is one such approach 

HSIiiEi- The 

Pauli potential is a non-physical repulsive potential in- 
troduced to avoid the overlap of the wave packets in the 
phase space during the time evolution. However, when 
it is applied in heavy ion collisions, the Pauli potential 
operates as a spurious repulsive force to increase the nu- 
clear stopping, especially during the early stages of the 
collisions Recently, Papa et al. introduced a new 

procedure to remove such overlaps of the wave packet in 
the phase space at each time step, without introducing 
the Pauli potential. The model is called the Constraint 
Molecular Dynamics model Although the idea is in- 
teresting and the computation time is short, no extensive 
application has yet been made for comparisons with ex- 
perimental results. 

In order to resolve the problem from the quantum me- 
chanical side, the Fermionic Molecular Dynamics model 
(FMD) and the Antisymmetrized Molecular Dynamics 
(AMD) model have been proposed 0, 0|. In both 
models the total wave function of the system is antisym- 
metrized and described by a Slater determinant of Gaus- 
sian wave packets. The time evolution of the centroid 
of the wave packets is treated in a classical manner. In 
FMD the width of the wave packets is treated as variable 
in time and NN collisions are treated as potential scat- 
terings. Until now the calculations have only been made 
with a harmonic oscillator potential and no application 
to heavy ion collisions has yet been made [l^ . 

In AMD-V, an improved version of AMD used in this 
study, the quantum nature of nuclcons in the wave packet 
propagation is incorporated as follows 0, 0| : 

• The total wave function of the system is antisym- 
metrized and therefore the Pauli principle is re- 
spected at all times. 

• The Pauli blocking in stochastic NN collisions is 



taken into account in an unambiguous manner. 

• The probabilistic nature of the wave packet is taken 
into account as a diffusion process during the wave 
packet propagation. The diffusion process is formu- 
lated in a manner to take into account the quantum 
branching to many final states of multifragmenta- 
tion channels. 



AMD-V has been applied for intermediate heavy ion 
collisions and found to reproduce reasonably well the ex- 
perimental results [13, UJI ISflj . For example in the pre- 
vious study of the ^'^Zn + ^^Ni reactions at 35A-79A 
MeV 10 , the calculated multiplicity, charge distribution 
and energy spectra of the reaction products were in good 
agreement with the experimental results. In that study 
it is pointed out that nuclear semi-transparency plays an 
important role in the multifragmentation process. In the 
present paper we extend the earlier study with both ex- 
periments and model calculations on additional systems. 
Experimental reaction measurements have been extended 
both to heavier systems and to lower incident energies. 
For the AMD-V calculations, the effective interaction and 
in-mcdium nucleon-nucleon (NN) cross section are the 
two important ingredients. In this work two different ef- 
fective interactions with different stiffnesses of the equa- 
tion of state (EOS) and two different formulations of the 
in-mcdium NN cross sections have been employed. 

A goal of the present work is to elucidate the reaction 
mechanisms, especially focusing on the equilibration and 
multifragmentation processes, using a reliable dynami- 
cal model for intermediate heavy ion reactions. For this 
purpose, the best parameter set for AMD-V calculations 
is chosen by comparing different calculations to the ex- 
perimental results. Then employing the best parameter 
set in AMD-V, the reaction mechanisms are studied in 
detail for the calculated events. Along these lines this 
paper is organized as follows; In Sec. II the experiment 
is described. In Sec. Ill a brief description of the AMD 
approach, effective interactions and in-medium NN cross 
sections is presented. In Sec. IV some remarks on the 
data analysis both for experimental and calculated re- 
sults are given. In Sec.V detailed comparisons between 
the experimental results and calculated results are pre- 
sented. In Sec. VI, a detailed analysis of the underline 
reaction mechanisms is presented and possible multifrag- 
mentation scenarios are proposed. In Sec. VII, the pro- 
posed multifragmentation mechanisms and existing data 
are discussed. In Sec. VIII a summary is given. 

In this paper nine reactions, ^^Zn + ■'^*Ni, ^^Zn + ^^Mo, 
^'^Zn + Au, at 26A, 35A and 47A MeV, have been 
studied. In order to refer to each reaction system in the 
text, the target name and the incident energy are used 
throughout the paper for simplicity, thus the ^'*Zn -I- ^*Ni 
reaction at 47A MeV becomes ^^Ni at 47A MeV. 
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TABLE I: NIMROD Charge Particle Array 


Ring 


Angle 


No. of segments 


Solid Angle 




(deg.) 




(msr) 


1 


4.3 


12 


0.96 


2 


6.4 


12 


2.67 


3 


9.4 


12 


4.26 


4 


12.9 


12 


7.99 


5 


18.2 


12 


16.1 


6 


24.5 


24 


12.7 


7 


32.1 


12 


33.6 


8 


40.4 


24 


27.6 


9 


61.2 


16 


154.0 


10 


90.0 


14 


207.0 


11 


120.0 


8 


378.0 


12 


152.5 


8 


241.0 




II. 


EXPERIMENT 





The experiment was performed at the K-500 super- 
conducting cyclotron facihty at Texas A&M University, 
using the An detector array, NIMROD, (Neutron Ion 
Muhi-detector for Reaction Oriented Dynamics). ^"'Zn 
projectiles were incident on ^^Ni, ^^Mo and ^^^Au tar- 
gets at energies of 26A, 35A and 47A MeV. NIMROD 
consists of a charged particle array set inside a neutron 
ball. The charged particle array is made of 166 segments 
in 12 concentric rings around the beam axis. Eight for- 
ward rings have the same geometrical design as the IN- 
DRA detector, but have less granularity '2l|. The angle, 
number of segments in each ring and solid angle of each 
Csl segment are given in Table I. 

The eight forward rings are covered by ionization 
chambers (IC). Furthermore in each of these rings two of 
the segments have two Si detectors between the IC and 
Csl detectors (super telescopes) and three have one Si 
detector. Each super telescope is further subdivided into 
two parts. The Csl detector is a Tl doped crystal read 
by a photo-multiplier tube. A pulse shape discrimination 
method is employed to identify particles, using different 
responses of fast and slow components of the light output 
of the Csl crystals for different charged particles . The 
ionization chambers were made of fiber-glass (GIO) and 
were filled with 30 Torr of CF4 gas. Front and back win- 
dows were made of 2.0 /xm aluminized Mylar foil. The 
signals were read by 5 to 10 fine wires, arranged per- 
pendicular to the particle direction in the active volume. 
In the Csl detector Hydrogen and Helium isotopes are 
clearly identified and Li fragments are also isolated from 
the heavier fragments. In the super telescopes, all iso- 
topes with atomic number Z < 8 are clearly identified 
and in all telescopes, particles are identified in atomic 
number. 

The energy calibration of the Si detectors was made 



with a ^^^Th source and the observed punch through en- 
ergies of identified particles. The punch through energies 
are calculated using a Range- Energy table . Since the 
energy losses of light particles, especially high energy Hy- 
drogen isotopes, are rather small in the Si detectors, eval- 
uation of the energy deposited in the Csl crystal requires 
special care. Therefore an additional energy calibration 
was performed as a separate run using a few telescopes 
in an 80 cm diameter scattering chamber. In the cali- 
bration run, the reaction of '^'*Zn -I- ^^Mo at 47A MeV 
was chosen as the standard. Si-Si telescopes backed by 
Csl detectors of three different lengths (1 cm, 3 cm and 
5 cm) were used to measure the inclusive energy spectra 
of light charged particles. The energy spectra were mea- 
sured at all angles corresponding to those of the 12 rings 
and solid angles were adjusted to be similar to those in 
Table I at each angle. The energy calibrations for high 
energy particles were made using the punch-through en- 
ergies of different lengths of the Csl crystals. The energy 
calibration for heavier fragments were made, using the 
Si detector calibration. The extracted energy spectra of 
fragments with atomic charge between 4 to 10 have been 
compared with those in reference for the ^''Zn + ^^Ni 
reaction at 35A MeV, in which Si detectors were used as 
AE-E telescopes, and a good agreement is obtained. 

Neutron multiplicity was measured with the Air neu- 
tron ball surrounding the charged particle array. The 
neutron ball consists of two hemispherical end caps and 
a central cylindrical section. The hemispheres are 150 
cm in diameter with beam pipe holes in the center and 
they are upstream and downstream of the charged par- 
ticle array. The central cylindrical section is 1.25m long 
with an inner hole of 60 cm diameter and 150 cm outer 
diameter. It is divided into 4 segments in the azimuthal 
angle direction. Between the hemispheres and the central 
section, there are 20 cm air gaps for cables and a duct 
for a pumping station. The neutron ball is filled with 
a pseudocumene based liquid scintillator mixed with 0.3 
weight percent of Gd salt (Gd 2-ethyl hexanoate) [23 |. 
Scintillation from a thermal neutron captured by Gd is 
detected by five 5-in phototubes in each hemisphere and 
three phototubes in each segment of the central section. 

In the experiment, data have been taken in two differ- 
ent trigger modes. One is the minimum bias trigger in 
which at least one of the Csl detectors fired. The other is 
the high multiplicity trigger which required that at least 
3-5 Csl detectors fired. The minimum bias trigger was 
scaled down, typically by a factor of 10, to reduce the 
rate of peripheral events. In order to reduce the neutron 
background, the beam was swept away in the injection 
line between the ECR source and the K-500 cyclotron for 
1 msec following detection of an event. 
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III. MODEL CALCULATIONS 



and Qij is defined as 



A brief description of the AMD model, including recent 
improvements related to the present work, is given in this 
section. Two important ingredients, effective interactions 
and in-medium nuclcon-nucleon (NN) cross sections, are 
also described. In order to show possible effects on the 
experimental observables resulting from different ingre- 
dients, changes in nuclear semi-transparency, a charac- 
teristic feature predicted by the model for intermediate 
heavy ion collisions, are explored. 



A. AMD-V model 

In AMD a reaction system with N nucleons is described 
by a wave function which is a single Slater determinant 
of N Gaussian wave packets , 



d 



^{Z) = det exp|-zy(^rj 



^Z?}x.,(j) 



, (1) 



where the complex variables Z = {Z^; i — 1, . . . ,N} — 
{Zia-; i = 1, . . . , TV, a = x,y, z} represent the centroids 
of the wave packets. can be described by 



K, 



(2) 



The width parameter v is taken as = 0.16 fm-2 and 
Xoii represents the spin and isospin states of p "f, p J,, 
ri "f, or n J,. For a dilute nuclear gas system, and 
Ki correspond to the position and momentum of each 
nucleon. Inside the nucleus, however, these quantities 
do not have physical meanings because of the antisym- 
metrization. The time evolution of Z is determined by 
the time-dependent variational principle and the two nu- 
cleon collision process. The equation of motion for Z 
derived from the time-dependent variational principle is 



dZjr an-L 



> . ^-.j- dt ~ dzi 



(3) 



CiajT is a hermitian matrix defined by 



dz^dz,^ 



log($(Z)|$(Z)), (4) 



and TL is the expectation value of the Hamiltonian after 
the subtraction of the spurious kinetic energy of the zero- 
point oscillation of the center-of-masses of fragments. 
Two nucleon collisions are introduced by the use of the 
physical coordinates W= j" which are defined as 



(5) 



5(2* • Z,) 



log($(Z)|<f(Z)). 



(6) 



In molecular dynamics models with Gaussian wave pack- 
ets, the i-th nucleon at time t=to is represented in phase 
space by 



/,(r,p,to) = 8exp|-2i^(r-R,(io))' 



(p-P.(to)) 



2h' 



with the centroid and P.^. The total one-body dis- 
tribution function is the sum of fi. In AMD, this repre- 
sentation of nucleon as a simple Gaussian wave packet is 
only approximately valid when the physical coordinate 



2%^ 



(8) 



is used for the centroids. 

In order to properly treat reactions with many branch- 
ing channels, such as multifragmentation processes, 
AMD has been extended by introducing the wave packet 
diffusion effect as a quantum branching process. This ex- 
tended AMD is called AMD-V, since the wave packet dif- 
fusion effect is calculated with the Vlasov equation [l^ . 
The AMD-V code has been further improved in order to 
save CPU time in the numerical calculations and to be 
applicable to heavier reaction systems. In the newly de- 
veloped code, used for all calculations in this paper, the 
calculation of the wave packet diffusion effect has been 
reformulated and a triple-loop approximation has been 
incorporated [l^ . 

The calculations were performed in the VPP700E su- 
percomputer facility in RIKEN, Japan. For the ^^Zn -f 
^^Ni and '^''Zn + ^'^Mo reactions, about 5000 and 3000 
events, respectively, were generated at each energy in the 
impact parameter range of 0-12 fm for a given parameter 
set. For the ^'*Zn -f ^^^Au case, about 1000 events were 
generated in the impact parameter range of 0-14 fm. The 
calculations were started with a distance of 15 fm along 
the beam direction between centers of the projectile and 
the target. The calculation for each event was carried 
out typically up to t=3GG fm/c for the reactions at 35 A 
and 47A MeV and up to 500 fm/c for those at 26 A MeV. 
In the following text, t=0 is set as the time at which the 
projectile and the target touch each other. 



B. Effective interaction and in-medium NN cross 
section 

The Gogny interaction |2g| has been used successfully 
in the previous analyses 0, 0, 0, H^l ■ This interac- 
tion gives a soft equation of state (EOS) with an incom- 
pressibility value K of 228 MeV for infinite nuclear mat- 
ter and has a momentum dependent mean field. In the 
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TABLE II: Parameters for the stiff effective interaction 



dk 


Wk 


Bu 


Hk 


Mk 


a tp 


^surf 


[fm] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV fm^(i+"'] 


[MeV fm^] 


fc = l 0.7 


-402.4 


-100.0 


-496.2 


-23.56 






fc = 2 1.2 


0.96 


51.575 


46.535 


-33.41 
















1.24 1896 


75 



present work, a modified Gogny interaction with a stiff 
EOS with K=360 MeV is also appHed Hlj. The parame- 
ter set adopted for K=360 MeV was that used by Haddad 
et al. and labeled as D1-G3. However this new in- 
teraction gives neither the proper charge radius nor the 
correct binding energy for the nuclear ground state in 
AMD. In order to get proper ground state properties of 
nuclei, it is necessary to add a surface term in the Hamil- 
tonian. Furthermore we correct the two body interaction 



term, so that, the force produces a reasonable equation 
of state for asymmetric nuclear matter. The expectation 
value of our stiff Gogny force V is given by 

(V) = C£v2{r^,r,)) + (^t,p(r„r,)) + Vs„.f, (9) 

i<j i<j 

where 



-(r,-r,)Va"fc 



fc=1.2 



Wp(r.i,rj) 

Vsurf 



tp{l + P,)p{r,yd{r,~r,), 
Uuri J d3r^(a/3|F<,|a/3)[Vp„(r)] • [Vp^(r)]. 



(10) 

(11) 
(12) 



a,l3 



r 



The indices a and P take four states of spin and isospin, 
p 'I , p I, n 'I , and n J,. The parameters are shown in 
Table UTI This new effective interaction is used in all cal- 
culations with the stiff EOS in the present study. The 
incompressibility K is still 360 MeV, because the surface 
term does not contribute to the incompressibility of in- 
finite nuclear matter and the correction to the two-body 
part does not affect the symmetric nuclear matter. 

Calculated binding energies for the soft EOS deviate 
from the experimental values, which were taken from [28j . 
by about 0.5 MeV/nucleon for the mass number A < 30 
and about 0.2 MeV/nucleon for the heavier fragments 
with A < 100. The agreement to the experimental values 
become better for the heavier fragments. The calculated 
binding energies for the stiff EOS, on the other hand, 
show larger deviations and they are up to 1 MeV/nuclcon 
less bound for all fragments, comparing to the experimen- 
tal values. 

In the previous analyses, an empirical in-medium 
nucleon-nucleon cross section was used 0, 0, . This 
cross section is given by [l^ 



100 



1 + E/{200MeV) + 2min((p/ poY^^ , 1) 



(13) 



where p is the nuclear density and Po is its normal value. 
The cross section is given in mbarn. The calculated 
cross section for the normal nuclear density is shown in 
Fig. n^a) by a dashed line. In the cross sections no dis- 
tinction was made between p-p(n-n) and n-p collisions. 
The density dependence is taken into account, but the 
dependence is rather small in the range of the density 
variation expected for the reactions studied in this work. 
For p/pq = 1.5, the cross section decreases by about 10%. 

Li and Machleidt calculated the in-medium NN cross 
sections based on the Dirac-Brueckner approach js^, |^ . 
For np collisions, the in-medium cross section is given by 



fl.53|2.9i 



apn = [31.5 + 0.092 x [20.2 - E' 

1.0-h0.0034£;i-5V^ 
^ 1.0-t-21.55pi-34 

and for pp, the cross section is given by 

app = [23.5 + 0.00256 x [18.2 - £;"-5|4.0] 

1.0-h 0.1667£;ioV 
^ 1.0-h9.704pi-2 



(14) 



(15) 



A normal nuclear density of 0.18 fm is assumed in the 
paper. The calculated cross sections for normal nuclear 
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matter are represented by symbols in Fig. ^a). At low 
energies these cross sections are about two to four times 
greater than those resulting from the empirical formula 
previously employed. Above 100 MeV, both pp and np 
cross sections become smaller than that of the empirical 
prescription. For p^l.bpo, the cross section decreases by 
about 35% at E=0 MeV and 25% at E=100 MeV. 

In Fig. n^b) the number of attempted and Pauli- 
allowed collisions are shown for actual calculations with 
the empirical cross section and with the Li-Machleidt 
cross section. The calculations were made for the central 
collision events ( b < 3 fm) of ^^Zn + ^^^Ni at 47A MeV. 
The soft Gogny interaction was used. For the empiri- 
cal formula, the number of attempted collisions reaches 
~8 collisions/ (fm/c) at a time of 30 fm/c when the two 
nuclei totally overlap and decreases quickly to 3-4 colli- 
sions/ ( fm/c). Only about 20% of the attempted colli- 
sions are Pauli-allowed during this process. For the Li- 
Machleidt cross sections, the number of attempted colli- 
sions increases by almost a factor of two. About 30% of 
the attempted collisions are Pauli-allowed in this case. 



C. Nuclear semi-transparency 

It has been reported that nuclear transparency plays 
an important role for the multifragrnentation process in 
intermediate heavy ion collisions [3j US ■ The stiffness 
of the effective interaction and the in-medium NN cross 
section are both important ingredients for determining 
the degree of the transparency in the calculations. In 
the present study three parameter sets have been inves- 
tigated for the calculations. They are 

• soft EOS + NNemp (empirical NN cross section) 

• stiff EOS + NNemp 

• stiff EOS + NNlm (Li-Machleidt cross section). 

In order to show how much nuclear semi-transparency 
changes with the different parameter sets, parallel veloc- 
ity distributions for all nucleons (free nucleons and nu- 
cleons bound in fragments) are shown in Fig.|21 at a time 
t = 280 fm/c, for central collisions in reactions at 47A 
MeV. In the case of ^^Ni (top row), the majority of the 
projectile and target nucleons move in the same direction 
after the collisions. For the soft EOS (left column) the 
projectile nucleons exhibit a broad distribution, centered 
at about half of the incident velocity after penetrating 
through the target nucleus, indicating significant trans- 
parency. The transparency is reduced for the stiff EOS 
as seen in the right column of the figure. It is interest- 
ing that only a small difference is observed between the 
results calculated with NNg^p (middle) and those with 
NNlm (right). Therefore, in the framework of the AMD- 
V, the transparency depends significantly on the stiffness 
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FIG. 1: (a) Calculated in-medium NN cross sections for normal 
density nuclear matter. The cross sections calculated by the em- 
pirical formula Eq. Iliji are depicted by a dashed line and those of 
the Li-Machleidt formulae Eqs. 1141 . 1151 are given by dots(np) and 
squares (nn,pp). (b) Number of collisions as a function of reaction 
time for central collisions of ^^Zn -|- ^^Ni at 47A MeV. Events with 
b<3 fm are analyzed. Solid symbols indicate the number of at- 
tempted collisions and open symbols indicate the number of Pauli- 
allowed collisions. Circles show the results of the Li-Machleidt for- 
mulae and squares show those of the empirical formula. 



of EOS, but depends only weakly on the in-medium NN 
cross sections. As the result the parallel velocity distri- 
butions of all nucleons show a two-peak structure for the 
soft EOS and these two peaks tend to merge into one 
peak for the stiff EOS. When the target becomes heav- 
ier, the two-peak structure for the soft EOS is also less 
prominent and becomes one peak for the ^^^Au target. 
In that case no notable difference is observed between 
the soft and stiff equation of state. The differences of 
the peak velocities between nucleons from the projectile 
and those from the target are also plotted in Fig. O En- 
hancement of the semi-transparency for the soft EOS is 
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clearly seen as the sharp increases of the differences for 
^^Ni and ^^Mo. 



D. Afterburner and switching time 

The fragments generated in AMD-V calculations are 
generally in an excited state at a time of 280 fm/c. AMD- 
V should properly treat the cooling of these fragments in 
a quantum statistical manner |33l |. However in order to 
cool the fragments down to the ground state, a great 
deal of CPU time is needed. Instead of continuing the 
AMD-V calculation for such a long time, we stopped the 
calculation at t = 280 fm/c (which corresponds to a re- 
alistic CPU time to get a few thousand events in the 
VPP700E) and the fragment cooling was followed, using 
a statistical decay code as an afterburner. A modified 
version of GEMINI [^l was used as the afterburner. In 
this modified version, discrete levels of the excited states 
of light fragments with Z < 15 are taken into account 
and the Hauser-Feshbach formalism is extended to the 
particle decay of these fragments when the excitation en- 
ergy is below 50 MeV. Each AMD-V event was used 100 
times in the afterburner in order to sample all possible 
decay paths of the excited fragments. This also provides 
enough statistics for detailed comparisons to the exper- 
imental results. The switching time of t = 280 fm/c 
was chosen only for practical reason of the computation 
time in the VPP700E. In the AMD-V calculations, later 
switching times are preferable, because the particle evap- 
oration occurs in the quantum statistical manner [s^ . As 
discussed in reference ^3 the switching time of t = 280 
fm/c is late enough so that the final results do not depend 
significantly on this choice. 
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FIG. 2: Calculated parallel velocity distributions of nucleons at 
280 fm/c are plotted in the center of mass system for ^*Ni, ''■^Mo 
and ^^^Au at 47A MeV from top to bottom, respectively. The 
results for soft EOS -I- NNemp, stiff + NNemp and stiff + NN/,m 
are plotted from left to right, respectively. Thick solid and dashed 
lines indicate the contributions of nucleons from the target and 
the projectile, respectively. Thin lines show the sum of these two 
contributions. 



IV. DATA Analysis 

In order to perform direct comparisons between the ex- 
perimental results and the calculations, the efficiency of 
the experimental acceptance, such as the neutron ball 
efficiency and multi hit effects, have to be evaluated. 
Event classification is also crucial for the comparisons, 
because many observables change drastically depending 
on the impact parameter. In this section these experi- 
mental conditions and event classification are described 
in detail. 
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A. Neutron ball efficiency 

Neutron balls have been applied to measure neu- 
tron multiplicity in heavy ion reactions for the last two 
decades In order to simulate the neutron ball ef- 

ficiency the program DENIS has been widely used [s^ . 
This code is designed for low energy neutrons and no 



FIG. 3: Difference of the peak values of the parallel velocity dis- 
tribution of nucleons from the projectile and from the target. Peak 
values are obtained by a Gaussian fit around each peak. Dashed 
lines are to guide the eye. 
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secondary neutron generation is taken into account. Re- 
cently Trzcinski et al. developed a more sophisticated 
program MSX which takes into account secondary neu- 
tron generation [s^. These codes have been applied to 
neutron energies up to a few tens of MeV for neutron 
balls containing a minimum amount of material inside 
the scattering chamber. However NIMROD has a large 
amount of material inside the chamber. This can result in 
scattering, absorption and generation of neutrons. To ad- 
dress these possibilities we have used the GCALOR code 
coupled to the GEANT-3 simulation package to simulate 
the neutron ball in NIMROD [Ullg. The GCALOR 
code is designed to simulate low energy neutrons. It 
makes use of the latest compilation of cross section data. 
In the program the geometry and material of the Neu- 
tron Ball and the charged particle array have been taken 
into account in detail. Only cables, phototubes and Si 
detectors are neglected. The calculations were made for 
neutrons emitted isotropically from the target with en- 
ergies from 2 to 100 MeV in 2 MeV steps. For each 
event, one neutron is emitted at the target and one mil- 
lion events were generated for each case. Each neutron 
was followed until either it escaped from the neutron ball 
or became thermalized ( E < 0.03 eV). If it thermalized 
inside the liquid scintillator, we assumed the neutron was 
captured by the Gadolinium. When additional neutrons 
are generated, all neutrons were followed until they were 
thermalized or escaped. The detection efficiency for a 
neutron of a given initial energy is determined by the 
ratio of all detected neutrons to primary neutrons. The 
light detection efficiency inside the neutron ball was also 
simulated separately by generating 7 rays emitting from 
a point source scattered inside the neutron ball. For a 
threshold of 300 KeV the detection efficiency of a cap- 
tured neutron was nearly 100%. The program was run 
in two modes, i.e., with the charged particle array and 
without the array to specify different contributions to the 
efficiency. 

In Fig. 01 the calculated results for the neutron ball effi- 
ciencies are shown. For comparison the DENIS code pre- 
dictions are also shown by a solid line for the case with- 
out the charged particle array Below 10 MeV both 
calculations agree reasonably well. For the higher en- 
ergy neutrons, the efficiency calculated by DENIS drops 
below the results for the first detected neutrons in the 
GEANT simulation (triangles), because there is no neu- 
tron generation in DENIS. In other words, even for the 
first captured neutrons in the GEANT simulation, gen- 
erated neutrons start to contribute to the detection effi- 
ciency above 10 MeV, as one can see in the figure. This 
indicates that a significant number of neutrons are gen- 
erated in the liquid scintillator for these high energy neu- 
trons. The contribution from the generation of more than 
two neutrons becomes significant above 20 MeV. Since 
the contribution of the third detected neutron is only 
10% of that of the second detected neutron (not shown), 
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FIG. 4: Calculated neutron ball efficiency without the charged 
particle array (a) and with the charged particle array (b). In- 
trinsic neutron efficiency at a given energy is shown by dots. 
The contribution of the first detected neutrons is shown by 
triangles and that of the second and higher order detected 
neutrons is shown by squares. The solid line in the left figure 
shows results from the DENIS code. The open circles in the 
right figure display the detection efficiency of generated neu- 
trons when a proton of an initial energy, given on x-axis, is 
emitted at the target. 



the main contribution comes from one or two neutron 
generation in this energy range. With the charged parti- 
cle array, the contribution of generated neutrons almost 
doubles. The total efficiency of the neutron ball slightly 
increases above 20 MeV as the neutron energy increases 
and becomes even approximately constant with respect 
to energy. 

Protons also generate secondary neutrons at high en- 
ergy. However the generation of neutrons by protons is 
much less efficient than that by neutrons. This is because 
protons lose their energy very quickly by ionization pro- 
cesses and the cross section for neutron generation de- 
creases rapidly as the proton energy decreases. Only a 
5% contribution to the neutron efficiency is observed at 
a proton energy of 50 MeV and 14% at 100 MeV. 

The efficiencies shown above are those averaged over 
all angles. However in the intermediate heavy ion reac- 
tions the angular distribution of the neutrons has a sig- 
nificant forward peak, especially for higher energy neu- 
trons. Since the charged particle array has more material 
at forward angles and two air gaps between the central 
part and the two hemispheres are not symmetric rela- 
tive to the target position, the angular dependence of 
the neutron efficiency has to be taken into account for 
the actual application. Therefore the neutron efficiency 
was calculated as a function of both neutron energy and 
polar angle. The azimuthal angular dependence was ne- 
glected. Since the proton multiplicity is much smaller 
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FIG. 5: (a) Typical experimental two dimensional plot of fast vs 
slow components of charge integrated light output from a Csl de- 
tector. Insert shows the expanded spectrum for Hydrogen isotopes. 
Corresponding elements are indicated in the figure. HF stands for 
heavy fragments with Z > 5. (b) Distributions of artificially gener- 
ated double-hit events for pp, pa and aa. Each particle combina- 
tion is indicated in the figure. Solid lines correspond to the locus 
of the ridges of p, d, t, i^He, a, ®Li and ^Li from (a). 
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FIG. 6: (Left) Number of particles generated by the calculation 
for ^^Mo at 47A MeV. Total number of charged particles is given 
by a histogram as a function of particle ID. The particle ID is 
given by ID = 0, 1, 2, 3, 4 for p, d, t, ^He, a and ID=Z+2 for Z>3, 
respectively. The number of single hit particles is given by dots and 
that of double hit particles is shown by squares. The number of the 
detected single hit particle is shown by crosses. (Right) Number of 
double a hit events is shown as a function of multiplicity for the 
experiment (dots) and the simulation (histogram). In both cases 
one million events were analyzed. The calculation is done for the 
stiff EOS -I- NN/^jv/ case. The calculated events have been treated 
by the experimental filter. 



than that of neutrons, the proton contribution to the 
neutron ball efficiency was not taken into account in the 
present analysis. 

B. Particle identification and multi-hit events 

In NIMROD, light charged particles with atomic num- 
ber Z < 3 are identified by a pulse shape discrimination 
method in the Csl detector. A typical two dimensional 
plot of fast versus slow charge integrated signals is shown 
in Fig. Ela). Each type of particle lies along a specific 
curve and one can clearly identify different particles as 
indicated in the figure. In the insert the spectrum is 
expanded for Hydrogen isotopes. The events along the 
far left line correspond to 7 rays (and accidental cosmic 
muons). Hydrogen isotopes are clearly identified down 
to a few MeV/nucleon. 

Special care has to be taken in the pulse discrimination 
method when more than one charged particle hits a Csl 
detector. In order to simulate double hit events in an 
actual detector, two experimental events were artificially 
mixed and the two outputs for a given combination of 
light charged particles were added both for the fast and 
slow components. The loci of double-hit events for pp, 
pa and aa are plotted in Fig. Isfb). The pp events are 
scattered along the deuteron and triton lines. Therefore 
these double hits are identified as a single deuteron or 
a triton in the present experimental data analysis. The 
main part of pa events lie near the alpha line. Most of 



these events, therefore, are identified as a single a par- 
ticle. When two a particles hit a detector, these events 
lie between the a and ^Li lines and are easily identified 
as two a particles, though only their summed energy is 
given experimentally. This is clearly seen in the experi- 
mental spectrum in Fig. Efa) and indicated by 2a. 

The rate of multiple hits has been simulated using 
AMD-V. Results are shown in Fig. Eta) for the case of 
^■^Mo at 47A MeV. The total number of each generated 
particle, the number of single hits and that of double hits 
are given by different histograms. For the double hits 
both of the particles are counted separately. According 
to this simulation, about 12% of protons and 20% of a 
particles hit a Csl detector which is hit by other charged 
particles. The rate of two a particles is larger, mainly be- 
cause of the two a decay of ^Be. The rate increases as the 
detector angle increases, because the solid angle of each 
segment increases rapidly as the angle increases as seen 
in Table I. The actual detected number of each particle 
as a single hit, after filtering through the experimental 
filter, is also given in the figure. About 85% of single-hit 
light charged particles were detected and identified. 

In order to verify the double-hit rate in the above sim- 
ulation, a comparison was made between two-a events in 
the experiment and in the simulation. The results are 
shown in Fig.E^b). A reasonable agreement between the 
experimental and calculated results is obtained. About 
10% of the events have one double a hit, both in the 
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FIG. 7: Calculated effects of double hits on the energy spectra of 
deuterons (left) and a (right) particles at three different detection 
angles, as indicated in each figure. Reaction for ^■^Mo at 47A MeV 
is used. Thick line histograms show the single hit spectra and thin 
line histograms indicate the spectra of pp in the deuteron spectra 
and pa in the a spectra. The energy of the double hit events is 
given by the sum of the energy of the two particles. 



experiment and in the simulation. 

The effect of muhi-hit events on the energy spectra of 
Hght charged particles has also been studied. In Fig. [7| 
deuteron and alpha spectra at three different detection 
angles are shown. In the calculation all pp-double hits 
are assigned as a deuteron hit with the summed energy 
and all pa hits are assigned as a single a hit. The con- 
tribution of the double hits at 9 = 9.4° is a few percent 
in both cases, and gradually increases as angle increases. 
At backward angles, especially for deuterons, the contri- 
bution becomes of the same order as that of the single 
particle hits. The contributions for the a spectra are 
slightly less. 

For heavier fragments with Z > 4, particle identifica- 
tion was made by the AE-E method and no multi-hit 
problem in particle identification occurs. (The energy 
loss of a proton or an a particle is much smaller than for 
a heavier fragment.) However heavier fragments from the 
target-like source have rather small kinetic energies and 
may be stopped in the AE detector. In this experiment 



heavy particles emitted at angles larger than 9 = 45° 
were not identified by charge. As seen in Fig.Efa) heavy 
fragments with Z > 3 were detected with an efficiency of 
40-50%. 

In the experimental filter applied to the calculated 
events, all pp hits are identified as a deuteron, all pa 
are identified as an a and aa events are treated as two-a 
hits. A heavy particle ( Z > 3) with light charged parti- 
cles in one detector is identified as a single heavy particle 
hit and the light charged particles are not counted in the 
charged particle multiplicity. Events with more than two 
hits in a single detector are treated as not identified for 
the largest fragment with Z < 2 and as a single hit for 
the largest with Z > 3, though the rate of these events is 
very small in the reactions studied in this paper. 



C. Event classification 

Detected events have been classified in four groups (Vi- 
olent, Semi- Violent, Semi-Peripheral and Peripheral), de- 
pending on the violence of collisions. This assignment is 
based upon the neutron and light charged particle mul- 
tiplicities and the total transverse energy of the light 
charged particles. In Fig. |H1 typical two dimensional plots 
of the normalized total multiplicity M^p/ Asystem versus 
the normalized total transverse energy are 
shown for ^^Mo at 47A MeV, both for the experimental 
and the calculated results. In both cases the experimen- 
tal inclusive distribution (top) shows a slightly broader 
distribution than the distribution calculated with AMD- 
V. For the multiplicity axis, this is mainly induced by 
the neutron multiplicity distribution as seen in the next 
section. The top 5% (3%) of the minimum bias events in 
the experiments were assigned to "Violent" collisions for 
58Ni and ^^Mo (^^^Au). The next 20% (10%) events were 
assigned as "Semi- Violent" collisions and the following 
20% (10%) of events were assigned as "Semi-Peripheral". 
The rest were assigned as "Peripheral" collisions. The 
same boundaries were applied to the calculated events. 
The resulting distributions corresponding to these cuts 
are compared in Fig. |H1 

The distributions of impact parameter for the four dif- 
ferent class of events were studied using AMD-V and are 
shown in Fig. |51 For the "Violent" and "Semi- Violent" 
classes the calculated distributions are very broad and 
about 50% of events of each class overlap with the distri- 
bution of another class. In the "Violent" class for ^®Ni, 
the distribution reaches up to 8 fm and about 70% of 
events are distributed in the impact parameter range of b 
< 5 fm. The distributions become broader for the heavier 
targets. In the "Violent" class for ^^^Au, the distribution 
reaches up to 10 fm and only about 35% of events origi- 
nate from collisions with b < 5 fm. On the other hand the 
events in the "Semi-Peripheral" and Peripheral classes 
show a rather localized distribution with a full width at 
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FIG. 8: 2D plots of the experimental (right) and the calculated 
(left) MLp/Asyste,n vs Ef C-P/Ej„,„ for ^^Mo at 47A MeV. M^p 
is the observed multiplicity of light particles, including neutrons, 
and Asystem is the total nucleon number of the reaction system. 
E^*-^^ is the sum of transverse energy of the light charged particles 
with Z < 2 and E^eam is the projectile incident energy. Gener- 
ated events by AMD-V has been filtered through the experimental 
acceptance. Contours are in logarithmic scale and the scale is set 
arbitrarily. 



a half maximum (FWHM) of 2-3 fm, indicating that most 
of these events actuahy originated from coUisions with a 
large impact parameter. Therefore throughout this pa- 
per, the word "Violent" is used instead of "Central" for 
the class of events with the highest multiplicity and the 
largest transverse energy. The word " Peripheral" is used 
as is customary. 



V. COMPARISON BETWEEN EXPERIMENT 
AND MODEL 

Detailed comparisons of the experimental and calcu- 
lated results are presented in this section for the direct 
observables, such as multiplicity distributions, charge dis- 



FIG. 9: Impact parameter distributions for different classes of 
events for (a) ^Sni (b) ^^Mo and (c) ^^^Au at 47A MeV. Trian- 
gles (down), triangles (up), squares and dots indicate the results 
for "Violent" , "Semi- Violent" , "Semi-Peripheral" and " Peripheral" 
classes, respectively. The area of each distribution is normalized to 
1. 



tributions, energy and velocity spectra. Comparisons 
for the elliptic flow are also shown. In order to make 
these comparisons, all calculated results have been fil- 
tered through the experimental conditions, unless other- 
wise specified. 



A. Multiplicity distributions 

Neutron and charged particle multiplicity distributions 
are presented in Figs. ^1 to E| In Fig. ^| experimen- 
tal neutron multiplicity distributions are compared with 
the calculated results for the "Violent" collisions of the 
reactions at 47A MeV. In each figure the experimental 
and the calculated results for the three different sets of 
parameters are shown. Overall the experimental results 
show broader distributions than those of the calculations 
for all reactions. For ^^^Au a shoulder is observed above 
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FIG. 10: Neutron multiplicity distributions for events in the "Vi- 
olent" class of events of the reactions at 47A MeV for (a) ^*Ni, 
(b) ^^Mo and (c) '^^''Au. Experimental results are shown by cir- 
cles and calculated results are shown by different lines. Thin solid, 
dashed and thick solid lines indicate the results of Soft-|-NNemp, 
Stiff -|-NNemp and Stiff -t-NN^Mi respectively. No efficiency and 
background corrections were applied for the experimental distri- 
butions, whereas the calculated results have been treated with the 
experimental filter. All distributions are normalized to one million 
events in total. 



the neutron multiplicity M„ > 40. It originates from 
pile-up in which two reactions occur during the lO/usec 
beam period before the shut-off of the beam. The pile-up 
is related to the reaction rate and is different for differ- 
ent reaction systems and incident energies (it can also 
fluctuate with time). Therefore in this work the pile-up 
effect is not taken into account in the experimental filter. 
The mean values of the calculated distributions for ^^Ni 
and ^^Mo agree with those of the experiments within 1-2 
neutrons, whereas that for ^^''Au produces ~3 additional 
neutrons for the soft EOS and ~5 additional neutrons for 
the stiff EOS. While the mean value for ^^Ni and ^^Mo 
is slightly better reproduced by the calculations with the 
stiff EOS, that for ^^^Au is better fit with the soft EOS. 
For all reactions the calculated multiplicities for the stiff 
EOS, but different NN cross sections, show almost iden- 
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FIG. 11: Multiplicity distributions of light charged particles (Z 
<2) (left) and heavier charged particles (Z >3) (right) for events in 
the "Violent", "Semi- Violent" and "Semi-Peripheral" classes, from 
top to bottom, respectively, for ^-^Mo at 47A MeV. Selected classes 
are indicated in each figure. Experimental results are shown by cir- 
cles and calculated results are shown by different lines. Thin solid, 
dashed and thick solid lines indicate the results of Soft-|-NNemp, 
Stiff-(-NN emp and Stiff -f-NN^^ jv/ 5 respectively. All calculated results 
have been treated with the experimental filter. All distributions are 
normalized to one million events in total. 



tical distributions. 

Similar trends are also observed for other reactions. In 
the top-left of Fig. 1151 the mean multiplicities of neu- 
trons, corrected for the efficiency, are summarized for 
"Violent" collisions for all reactions studied here. The 
calculated results for Soft-t-NNemp and Stiff-|-NNiM are 
compared to the experiment. The experimental mean 
values for ^^Ni and ^^Mo are reasonably well reproduced 
by the calculations but for ^^^Au the experimental mean 
values are exceeded by about 10-20% at all incident en- 
ergies. 

Typical charged particle multiplicity distributions for 
light charged particles (Z < 2) and heavier fragments are 
shown separately in Fig. For collisions in "Violent" , 
"Semi- Violent" and "Semi-peripheral" classes for ^^Mo 
at 47A MeV, both experimental and calculated results 
are shown. Calculations with three different parameter 
sets lead to very similar distributions for all cases. For the 
"Violent" and "Semi- Violent" collisions the light charged 
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FIG. 12: Multiplicity distributions of selected particles for the 
"Violent" collisions for '■'^Mo at 47A MeV. Particles are indicated 
in each figure. See also the caption of Fig. 1111 



particle multiplicity is overpredicted by 20% in the cal- 
culations. The heavier charged particle multiplicity is 
also overpredicted by about one unit in the mean value. 
The widths of the distributions are well reproduced for 
all cases. For the "Semi-Peripheral" collisions, the light 
charged particle multiplicities are well reproduced both 
in the mean value and width, but for the heavier charged 
particle multiplicities all calculations show a shoulder at 
the multiplicities around 4-5. 

Typical individual charged particle multiplicities for 
the "Violent" collisions are shown in Fig.^jfor the same 
reaction. There is an interesting connection between the 
multiplicity of protons and a particles. For the soft EOS, 
the proton multiplicity distribution is reproduced quite 
well, both in the mean value and width. On the other 
hand the mean value of the a distribution is overpre- 
dicted by more than 2 units. On the contrary, for the 
stiff EOS the a multiplicity distribution is well repro- 
duced, but the proton multiplicity distribution is overpre- 
dicted by more than two units. Since the differences are 
small for the other light charged particles and their ex- 



FIG. 13: Proton multiplicity distributions for the "Violent" colli- 
sions of all reactions. Figures in the same column show the results 
from the same reaction system and figures in the same row show 
those at the same incident energy. The reaction system is indi- 
cated at the top of the figure and the incident energy is indicated 
in each row. Experimentally observed multiplicities are shown by 
the circles and calculations with soft EOS + NNemp and stiff EOS 
-I- NNij\/ are shown by thin and thick solid lines, respectively. 



perimental distributions are reasonably reproduced, the 
overprediction of the calculated multiplicity distributions 
for light charged particles seen in Fig. llll is mainly caused 
by alpha particles for the soft EOS and by protons for 
the stiff EOS. For Li fragments the calculated distribu- 
tions show slightly smaller values than the experimental 
values for all three calculations. Those for Be fragments, 
not shown, also show a similar trend. For heavier IMF 
with 5 < Z < 15, all calculations overpredict the experi- 
mental multiplicities by 1-2 units, as seen in Carbon and 
Neon cases in Fig.^] Here also one can see that the cal- 
culated results for the stiff EOS with different NN cross 
sections show almost identical distributions for all cases. 

Proton and a multiplicities for all reactions are shown 
in Fiefs. 1131 and 1141 for "Violent" collisions, respectively. 
For most reactions, calculations with the soft EOS are fa- 
vored for the proton multiplicity distributions and those 
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FIG. 14: Similar plot to Fig. Iliji but for a particles. 



with the stiff EOS are favored for the a multiphcity dis- 
tributions. For ^'^''Au at 26A and 35A MeV, however, 
the calculated proton multiplicities are overpredicted for 
both EOS's while the experimental a multiplicity distri- 
butions are rather well reproduced by the calculations 
with the soft EOS. 

For the "Violent" class the results are also shown in 
Fig. 1151 In the figure the efficiency corrected experimen- 
tal multiplicities and non-filtered calculated multiplicities 
are shown. The experimental detection efficiency for each 
particle is evaluated from the calculated events. The dif- 
ferences between calculations for a given particle agree 
within a few %. The error bars in the figure include the 
systematic errors. The increase of the experimental mul- 
tiplicity with energy and target mass is reasonably repro- 
duced for all particles by both calculations, whereas the 
calculated absolute mean values deviate from the exper- 
imental multiplicities by about 20-30% for some cases. 
For Li, the experimental mean values are overpredicted 
for most cases and the deviation becomes significant for 
"'^^''Au. For the heavier fragments the calculated mean 
values become larger than the experimental ones in gen- 
eral. For Ne, the experimental mean values are well re- 
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FIG. 15: Summary of mean multiplicities of selected particles for 
the "Violent" class of events for all reactions. Experimental values 
are efficiency corrected but calculated ones are non- filtered values. 
Reaction systems and incident energies are indicated on the x axis 
and particle is indicated in each panel. Experimental results are 
shown by solid dots and calculated results with soft EOS -I- NNemp 
and stiff EOS + NN^jv/ are shown by open squares and triangles, 
respectively. 



produced at 26A MeV, but overpredicted at 47A MeV 
by both calculations. 



B. Charge distributions 

Charge distributions for different classes of centrality 
for ^^Mo at 47A MeV are shown in Fig. [THl The exper- 
imental charge distributions evolve to the larger Z side 
when collisions become less violent. This trend is well re- 
produced by all calculated results, as one can see in the 
figure. No significant difference is observed between the 
calculated results with different parameter sets, except 
for fragments with 20 < Z < 30 in the "Semi-Peripheral" 
class. For most of the cases the calculated fragment mul- 
tiplicities with 5 < Z < 15 are overpredicted by a factor 
of 1.5-2. A similar discrepancy is observed at lower in- 
cident energies. In Fig. 1171 the charge distributions for 
^*Ni and ^^^Au are compared with the calculations. For 
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^^Ni the discrepancy is similar to ^^Mo and, for -'^^''Au, 
the discrepancy becomes larger. For ^^^Au all calculated 
results are overpredicted by a factor of 2-3. As one can 
recognize in the figure, the charge distributions for the 
two targets are quite similar both in the experiments and 
calculations. This is also true, in some extent, for all dif- 
ferent reactions studied in the present work. In Fig. 1181 
experimental efhciency corrected charge distributions are 
compared with each other for light fragments with 3 < 
Z < 20. The detection efficiency correction is made by 
comparing the calculated results with and without the 
experimental filter. Each group shows results for the 
three different targets at a given incident energy. For 
^^^Au fragment multiplicity tends to be slightly smaller 
at 35 and 47A MeV, but the shapes of the global dis- 
tributions are rather similar to each other, even for dif- 
ferent incident energies. A similar trend is also seen in 
the calculated events. The resemblance of the relative 
charge distributions is quite interesting because the re- 
action mechanisms are quite different as discussed later. 



C. Energy spectra 

For the violent collisions for ^^Mo at three different 
incident energies, typical energy spectra of light charged 
particles are shown in Figs.^jto Experimental spec- 
tra and calculated spectra for the soft EOS -I- NNg^p 
and the stiff EOS -I- NNlm are shown. Energy is scaled 
by the beam energy. The vertical axis is the differen- 
tial multiplicity and all results are given in an absolute 
scale. As one can see for all cases, the experimental spec- 
tra, at three incident energies, are very similar to each 
other in shape and angular dependence. This indicates 
that these energy spectra can be parametrized as emis- 
sion from moving sources with similar source velocities 
and apparent temperatures, but scaled by the incident 
energy. In the present work, however, no such analysis 
was performed. Instead the experimental energy spectra 
and angular distributions are compared with calculated 
AMD-V events. As seen in Fiars. 1191 to 1221 a reasonable 
agreement is obtained for all cases. For protons, however, 
the experimental slopes of the high energy tails in Fig. 1191 
tend to be harder than those in calculations, especially at 
26A and 35A MeV. The deviation becomes significant at 
angles of 6* = 40.4° and 61.2° which correspond to emis- 
sion near = 90° in the center of mass system. A similar 
discrepancy is also observed at 47A MeV. For the calcula- 
tions with the stiff EOS, an excess of low energy protons 
is clearly observed at angles between 24.5° < 9 < 60.2° 
at both 35 and 47A MeV. This is slightly less prominent 
at 26A MeV. This is the main cause of the excess of the 
proton multiplicity in the calculations with the stiff EOS, 
observed in Fig. ^1 

The experimental a energy spectra are compared with 
calculated spectra in Fig. 1201 At 35 and 47A MeV the 
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FIG. 16: Angle integrated charge distributions of events in (a) 
"Violent", (b) "Semi- Violent" and (c) "Semi-Peripheral" classes 
arc shown for ^-^Mo at 47A MeV. Experimental results are shown 
by circles. Dot-dashed, dashed and solid lines correspond to soft 
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experimental spectra are well reproduced by the calcu- 
lations with the stiff EOS, except at the most forward 
angle, whereas the calculations with the soft EOS over- 
predict the yields on higher energy side at most of the 
forward angles. This overproduction causes the overpre- 
diction of the alpha multiplicity observed for the soft EOS 
in Fig.[l3| On the other hand, at 26A MeV, the calcula- 
tion with the soft EOS reproduces the experimental spec- 
tra quite well at most angles. The stiff EOS significantly 
underestimates the yields in higher energy side. 

Deuteron and triton energy spectra are shown in 
Fig. Ell and Fig. [23 At 26A MeV both calculations un- 
derpredict the yields for both deuterons and tritons of 
higher energies. This is similar to observations for the 
proton spectra. For the deuteron spectra at higher inci- 
dent energies, the calculations with the stiff EOS repro- 
duce the experimental spectra better than those with the 
soft EOS at the four most forward angles. The overpre- 
diction of low energy deuterons, that is seen in proton 
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spectra, is also observed at angles between 24.5° < < 
60.2°. This excess is not observed for the triton spectra, 
but the corresponding statistical fluctuations are signifi- 
cantly large in that case. 

Typical inclusive energy spectra of heavier fragments 
are shown for Oxygen fragments in Fig. |23for ^^Mo at 
35 A and 47A MeV. As discussed in ref [l3j, the energy 
spectra of IMF depend little on the multiplicity class se- 
lection. The experimental spectra, scaled by the beam 
energy, show very similar shapes and angular distribu- 
tions for both incident energies, except at the most for- 
ward angle. This again indicates that these spectra can 
be described by moving source (or sources) parametriza- 
tions with similar source velocities and apparent temper- 
atures, scaled by the incident energy. In AMD-V most of 
experimental spectra are better reproduced by the calcu- 
lations with the stiff EOS, especially in the higher energy 
side. 



FIG. 17: Angle integrated charge distributions of events in violent 
collisions for (a) ^^Ni and (b) ^^''Au at 47A MeV. Dot-dashed 
and solid lines correspond to soft EOS -|- NNg^p and stiff EOS 
+ NNij\/, respectively. 
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FIG. 18: EfHciency corrected experimental charge distributions for 
fragments with 3 < Z <20. Each group shows results with three 
different targets at a given incident energy. The incident energy 
is shown on the right. The differential multiplicity is given in an 
absolute scale at 47A MeV and premultiplied by factors of 100 and 
10, respectively, for results at 26A and 35A MeV. 



D. Velocity distributions 

As shown in Fig. 13 for the parallel velocity distribu- 
tions of all nucleons at early stages of reactions, AMD-V 
predicts a distinct difference between the soft EOS and 
stiff EOS calculations for reactions with the lighter tar- 
gets at 47A MeV. The calculated distributions of all free 
and bound protons also show similar differences as seen 
in Fig. 1241 However, in order to make the distribution 
for all protons experimentally, all fragments have to be 
identified at all angles with very low energy threshold. 
This is a very difficult task in experiments. Instead we 
compare the experimental and calculated distributions of 
free protons. 

We first compare the calculated parallel velocity dis- 
tribution of all free and bound protons before the after- 
burner and that of only free protons after the afterburner. 
The results are shown in Fig.|^for the reactions at 47A 
MeV. The results with the soft EOS are shown on the 
left and those with the stiff EOS are shown on the right. 
For free protons after the afterburner the parallel veloc- 
ity distributions for the soft EOS show slightly broader 
distributions than those for the stiff EOS, but the shape 
becomes very similar. Especially for ^^Ni, on the left, 
the two-peak structure seen for free and bound protons 
before the afterburner becomes a single broad peak. As 
a result, no distinct difference is observed between cal- 
culated results with the soft EOS and the stiff EOS for 
the distributions of free protons after the afterburner . 
This indicates that the statistical decay process in the 
afterburner changes the distributions quite drastically. 

In the left column of Fig.|55lthe experimental parallel 
velocity distributions of free protons are compared with 
those of the calculations. The dips in the negative veloc- 
ity side in the experimental distributions are caused by 
the target shadow to the detectors at 90°. As one can see. 
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FIG. 19: Proton energy spectra for the violent collisions for ^■^Mo at three incident energies. The incident energy is indicated at the top 
of each figure. The spectra in each column correspond to the energy spectra at different angles, indicated in the left column. Experimental 
results are shown by dots. Thick and thin solid lines correspond to calculated results for Soft+NNe,7ip and Stiff+NN^jv/ , respectively. For 
all cases the differential multiplicity is given in absolute units. In order to avoid the overlap of the data spectra are multiplied by factors 
of 10" (n=0, 1,2, 3, 4, 5, 6, 7, 8, 9, 10) from bottom to top. No experimental filter is used for the calculated events. 



both calculations with the soft EOS and the stiff EOS re- 
produce the experimental proton distributions quite well 
and no distinct difference is observed between the two 
calculations. On the other hand the parallel velocity dis- 
tributions of a particles do show some differences be- 
tween the two calculations seen in the right column of 
Fig.l25lwhere the experimental results for a particles are 
compared with those of the calculations. In the calcu- 
lated results for ^^Ni, the two-peak structure is clearly 
observed in the simulation with the soft EOS, whereas 
a single broad peak is seen in that with the stiff EOS. 
This result is consistent with the calculations reported by 
Ono for ^OCa + *°Ca. at 35A MeV. This trend is also 
observed for ^^Zn -I- ^^Mo in Fig. I^S] but is slightly less 
prominent. For ^"^Zn -I- ^^^Au no difference is observed 
between the two calculations. The overall experimental 
trend for a particles favors the stiff EOS. 

Since the differences observed in the calculations with 
the soft EOS and stiff EOS in Fig. [3 are not reflected in 
the free proton distributions but appear in the a distri- 



butions, one can expect that the difference is related to 
the velocity distribution of heavier fragments. In Fig. 1251 
bi dimensional plots of parallel velocity versus perpendic- 
ular velocity are shown for fragments with Z=7,8 for the 
reactions at 47A MeV. The velocity distributions for cal- 
culations with the soft EOS show a clear two-peak struc- 
ture for both ^^Ni and ^^Mo. On the other hand those 
of the calculations with the stiff EOS show an elongated 
distribution. Although in the experiment about 50% of 
fragments are not detected because of the energy thresh- 
old and the detection angles, as shown by dashed lines in 
the calculated results, one can still clearly see that the ex- 
perimental IMF distributions favor the calculations with 
the stiff EOS. 



E. Collective flow analysis 

In collective flow analyses, the directed flow, az- 
imuthal angle correlation and elliptic flow have been 
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FIG. 20: Similar plots to Fig. 1191 but for a particles. 



studied [s^ ji^, 0| ■ These flows are essentially governed 
by particle emission at early stages. Therefore, the event 
multiplication in the afterburner does not reduce the sta- 
tistical errors on the results. Flow is also a rather small 
effect and a large number of events has to be generated 
in AMD-V to get reasonable comparisons to the experi- 
mental results. Of the three different flow analyses, the 
differential elliptic flow is an averaged observable and re- 
quires a smaller number of events for reasonable compar- 
isons. Therefore, we studied the differential elliptic flow 
V2 of protons ^4^], given by the averaged value of V2, and 
described by 



< V2{pt) >-- 



1 



N{pt 



(16) 



where N(pt) is the number of protons in a given pt range 
and Pjt is the transverse momentum of the j-th proton. 
Pjx and pjy are its perpendicular momenta in and out 
of the reaction plane, respectively. The reaction plane is 
determined, using transverse momentum of emitted par- 
ticles as proposed by Danielewicz and Odyniec In 
Fig. < V2 > values of the experimental and calcu- 



lated results are shown for protons for ^^Ni at 47A MeV. 
Events in the "Violent" and "Semi- Violent" classes are 
summed in order to increase statistics. Only protons in 
the mid-rapidity range (0.25 < y/jbeam < 0.75) are used. 
Errors are evaluated from the deviations from zero for 
< V2 > values when the azimuthal angle of the reaction 
plane is randomized. As one can observe, the experimen- 
tal < V2 > value shows a small deviation from zero in the 
energy range studied in this paper. No significant differ- 
ence is observed between the values extracted from the 
calculated events for different parameter sets. The cal- 
culated values are consistent with the experimental ones 
within error bars, except the most highest pt, where the 
statistics becomes poor. 



F. Discussion 

As seen in Fig. 12 a prominent difference is observed be- 
tween the calculated results with the soft EOS and with 
the stiff EOS in the parallel velocity distribution of all 
free and bound nucleons before the afterburner (t=280 
fm/c) for 5^Ni and ^^Mo at 47A MeV. The distribution of 
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FIG. 21: Similar plots to Fig. 1191 but for deuterons. 



all protons also shows a two-peak (or a shoulder) struc- 
ture as seen in Fig. |^ However this two-peak structure 
becomes a broad single-peak distribution for the free pro- 
tons after the afterburner. The experimental distribu- 
tions are well reproduced by both calculations with the 
soft EOS and with the stiff EOS. These good agreements 
suggest that the prominent difference, seen in the veloc- 
ity distribution at the early stages of the reaction, may 
be reflected in the bound nucleons in clusters. This is 
interesting because, in this case, the difference which ap- 
pears in the velocity distribution of the clusters may not 
depend on the details of the multifragmentation process. 
In other words, the signal of the stiffness of the EOS 
remains as a fossil signal in the velocity distribution of 
fragments. 

As seen in section V.B, our experimental multiplicity of 
fragments with Z > 4 is overpredicted in all calculations. 
A similar result has also been observed in the ^^®Xe-|-Sn 
reaction at 50A MeV . In order to resolve this discrep- 
ancy between the experimental and calculated results, an 
extension was made in the reference, so that the shrink- 
ing as well as the diffusion of the wave packets would be 
taken into account during their propagation in the mean 



field. This treatment with the new parameter suppresses 
the dynamical decay of the excited fragments at early 
stages and reproduces the experimental charge distribu- 
tion quite well. However the above observation for the 
parallel velocity distribution indicates that the signature 
of the stiffness of EOS remains as a footprint in the frag- 
ment velocity distributions with or without taking into 
account this kind of fine tuning of the model. 

In the present work, the Gogny interaction with the 
stiff EOS is favored by the experimental velocity distribu- 
tions of the clusters. On the other hand the Gogny force 
with the soft EOS gave a better description for the giant 
monopole resonance studies 01- However one should 
note the following: (1) In the giant resonance studies, 
normal nuclear matter is studied, whereas in this work 
the nature of compressed nuclear matter at high temper- 
ature is treated; (2) Our result strongly depends not only 
on the stiffness of the EOS but also on the reaction dy- 
namics, especially the Pauli blocking during the reaction. 
In AMD-V Pauli blocking is taken into account exactly 
for the propagation of the wave packets at all times. For 
stochastic nucleon-nucleon collisions, the Pauli-blocking 
is automatically taken into account in the transformation 
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FIG. 22: Similar plots to Fig. 1191 but for tritons. 



between the physical coordinate space and AMD space. 
However this transformation, which is given in Eq. jSJ, 
is an approximation. One of the consequences of this ap- 
proximation is seen in the charged particle multiplicity 
distributions for collisions near the surface, discussed in 
detail in ref ■ When the Pauli- blocking becomes very 
significant as in the reactions studied in this work, this 
procedure may cause some errors. In fact at 26A MeV the 
calculated results for both soft and stiff EOS significantly 
underpredict the slopes of energy spectra for protons, 
deuterons and tritons, as seen in Figs. IIHI2II22I This 
suggests that the AMD-V calculations lead to too much 
transparency at lower energies, resulting in less compres- 
sion and/or excitation energy at the early stage of the 
collisions. Therefore it is appropriate to only conclude 
that, in the framework of the present AMD-V calcula- 
tions, the stiff EOS is preferred to reproduce the present 
experimental results. 

Another interesting observation in this study is that 
the different NN collision cross sections do not alter the 
nuclear semi-transparency. Two formulations have been 
employed in the present calculations. Although the cross 
sections are different by a factor of 2-3 on average in the 



range of proton energies relevant to this study, the ap- 
parent effect on the nuclear stopping is rather small, as 
seen in Fig. [5] This result is quite surprising because in 
QMD studies of the nuclear collective flow, the effect of 
different NN cross sections on the balance energy is of 
a similar order to the effect of different stiffnesses of the 
EOS 113 . Ono et al. studied collective flow in the ''°Ar 
-I- ^^Al reaction, using AMD without the diffusion pro- 
cess In that study the strength of flow was reduced 
by about 30% at 45 A MeV when a, the empirical NN 
cross section given by Eq. ((T^ . is increased by 50%. 



The results of the differential elliptic flow analysis are 
consistent with the observation of the parallel velocity 
distributions, although the deviations of the experimen- 
tal and calculated results from the isotropic distribution 
are rather small. These results suggest that the diffusion 
process in AMD-V plays a significant role, not only in the 
multifragmentation process, but also in the wave packet 
propagation. This diffusion process tends to randomize 
the trajectory of wave packets and to smear out collective 
flow in this energy region. 
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FIG. 23: Inclusive energy spectra of IMF with Z=8 at different 
angles for violent collisions with ^■^Mo (a) at 35 A and (b) at 47A 
MeV. Experimental spectra are shown by symbols and calculated 
results with soft EOS + NNemp and stiff EOS + ^^lm are shown 
by thick and thin solid lines, respectively. The experimental spectra 
are plotted in an absolute scale, whereas each calculated result is 
normalized to the experimental spectra at S = 4.3°. 



VI. MULTIFRAGMENTATION MECHANISM 



FIG. 24: Calculated parallel velocity distributions of protons in 
the center of mass system for different reactions at 47A MeV. Each 
row corresponds to the same reaction which is specified by target 
as indicated in each panel on the left. The calculated results with 
soft EOS + NN emp are plotted on the left and those with stiff EOS 
+ NN^jv/ are on the right. The distributions of all free and bound 
protons before the afterburner at 280 fm/c are displayed by thin 
lines. The distributionis of free protons after the afterburner are 
shown by thick lines. The normalization for the latter is arbitrary. 



In the previous work on the ^''Zn + ^^Ni reactions at 
35A-79A MeV, we reported that the semi-transparency 
plays an important role for the multifragmentation pro- 
cess However for the ^"^Zn -t- ^^^Au reactions the cal- 
culations indicate that all projectile nucleons essentially 
stop in the target nuclei at all incident energies, as seen in 
Fig. 121 Therefore it is expected that transparency plays 
little role in the multifragmentation process for these re- 
actions and that other mechanisms, such as expansion 
and statistical multifragmentation processes, may play a 
dominant role in the disintegration of the system. On the 
other hand the similarity of the charge distributions for 
all reactions studied here, as seen in Fig.^| suggests that 
there is a common feature for the multifragment produc- 
tion in these reactions in which quite different dynamics 



are involved. 

Since the overall experimental results are reproduced 
reasonably well by AMD-V with the stiff EOS, in the fol- 
lowing we will analyze central collision events ( b < 3 fm 
) for the different systems, calculated with the stiff EOS, 
in order to elucidate the multifragmentation mechanism 
in great detail. No afterburner is applied in the analysis 
in this section. 



A. Global character of reactions 

In Fig. I2HI density distributions for ^^^Au arc plot- 
ted as a function of time. The projectile and the target 
are fully overlapped around t ~ 30 fm/c and the system 
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FIG. 25: Parallel velocity distributions of protons (left) and a 
particles (right) in the center of mass system for different reactions 
at 47A MeV. Reaction is specified by target as indicated in each 
panel on the left. Experimental results are shown by dots and 
the calculated results with the soft EOS +NNemp are plotted by 
histograms of thick lines and those for the stiff EOS + NN^m 
by histograms of thin lines. The calculated results are normalized 
to the experimental distributions to obtain the same number of 
protons at V// > 0. in each cases. 



starts to expand and undergoes multifragmentation for 
all incident energies. One should note that prefragments 
are already recognized at early stages of the reaction ( t 
~ 80^180 fm/c ). This is quite different from a statis- 
tical multifragmentation picture, in which a hot system 
expands and clusterizes at a low freezeout density |46l |. 
It should be also noted that the fragment sizes are very 
similar with each other for three incident energies. Nu- 
cleon emissions are identified as early as t 100 fm/c at 
35 A and 47A MeV and t - 200 fm/c at 26A MeV. 

In Fig. I^ni the mass number, density and excitation 
energy of the largest fragment (Fmax) and of the second 
largest fragment (F2, Z > 2) are plotted as a function 
of the reaction time for ^^^Au at 35A MeV. In order to 
evaluate these quantities, each calculated event is clus- 
terized in coordinate space with a coalescence radius of 
5 fm at a given time. In the first row, the mass numbers 
of Fjnax and F2 are shown. At t 200-300 fm/c, Fmax 
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FIG. 26: Two dimensional Galilean invariant plots of parallel ve- 
locity vs perpendicular velocity of all fragments with Z=7 and 8 for 
the reactions at 47A MeV. Plots are in the center of mass system. 
Each row shows the results for the same reaction system and reac- 
tion is specified by the target, indicated in the left figure of each 
row. The experimental distributions are shown on the left column. 
The calculated distributions for soft EOS + NNemp are shown in 
the middle column and those for stiff EOS -I- NN^jv/ are plotted on 
the right column. No experimental filter is applied for the calcula- 
tions, except the forward angle detection limit (3.5°). Dashed lines 
in the middle and right columns indicate the detection limit of the 
polar angle. Dashed circles indicate approximate energy thresholds 
for these fragment. 
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FIG. 27: Averaged differential elliptic fiow for protons for ^^Ni at 
47A MeV as a function of the transverse momentum pt- Events 
in the "Violent" and Semi- Violent classes are used. Experimental 
results are shown by dots and calculated results for stiff EOS -|- 
NN LM, stiff EOS -I- NNemp and soft EOS -I- NNemp are shown by 
open circles, squares and triangles, respectively. < V2 > values are 
calculated in 50 MeV/c steps in pt- Calculated results are plotted 
around the center value of each pt with a small shift to avoid overlap 
of results. 



undergoes multifragmentation. The generated F2 has a 
similar size to the Fmax after this multifragmentation. 

In the second row the maximum density Pmax is plot- 
ted. In the time of overlap, Pmax reaches around 1.4po 
and quickly returns to the normal density at t ~ 100 
fm/c. After that Pmax stays around 0.85po, Pmax of F2 
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FIG. 28; Time evolution of nuclear density distributions, projected 
on the reaction plane (x-z plane, z being the beam direction), for a 
calculated event with b ~ 2 fm at different reaction times for -'^^^Au. 
Incident energy is indicated at the top of each column. Reaction 
time is indicated in the unit of fm/c in each panel. The time zero 
is set at the time when the projectile and the target touch each 
other. The z axis is taken as the beam direction and the contour 
scale is in linear. The smallest circle indicates a nucleon. 



also shows a similar value to that of Fmax- This indi- 
cates that prefragments, as seen in Fig. |2H1 have density 
close to the normal density when they are formed at early 
stages and keep this density during the expansion of the 
system. 

In the third row, the density calculated from the root- 
mean-square radius (Rrms) is plotted. pRrms is cal- 
culated by assuming that the fragment has a spherical 
shape, i.e., 



PRrms Po 



^ r ^° 1 3 



(17) 



Here Aq, po and Ro are the mass number, density and 
r.m.s. radius of the initial projectile nucleus, respectively. 
PRrms reaches l.ipo at t ~ 30 fm/c , a similar value to 
Pmax, indicating that the shape of the composite system 
at the time of overlap is nearly spherical. pRrms decreases 
rapidly after that and reaches 0.2po at t ^ 200 fm/c, al- 
though pmax is close to thc normal density. This indicates 
that the system at this time has a very deformed shape 
and non-uniform density. pRrms of F2 is also small, indi- 
cating that emitted fragments have also an odd shape. 

In the fourth row the excitation energies are shown. 
The excitation energy of a fragmenti, E*, is calculated 
by subtracting the calculated binding energy from the 
internal energy, i.e., 



E* = Ej, 



E 



B.E. 



(T) + (V) + (Vc) - E 



B.E. 



(18) 



(T) is the expectation value of kinetic energy in the rest 
frame of the fragment, (V) is that of the effective inter- 
action given in Eq. 1^1 for the case of the stiff EOS, for 
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FIG. 29: Time evolution of mass number, density and excitation 
energy of the largest fragment (left) and the second largest fragment 
with Z I 2 (right) from top to bottom, respectively for ^^^Au at 
35A MeV. See details in the text. Contours are in logarithmic scale. 



example, (Vc) is that of the Coulomb interaction and 
Eb.b. is the calculated binding energy. The excitation 
energy of F^ax is E*/A ^ 5 MeV at the time of overlap 
and gradually decreases with time. An interesting ob- 
servation is that the excitation energy of F2 distributes 
around 3 MeV/nucleon independently on the emission 
time and shows significantly lower than that of F^ax 
at t < 300 fm/c. At t > 400 fm/c, the excitation en- 
ergy of F2 becomes similar to that of the F^ax- Some 
of F2 at t < 300 fm/c have excitation energies of less 
than 1 MeV/nucleon. These are Li isotopes, which have 
\Eb.e.\ < 5 MeV/nucleon. A more detailed discussion 
about cold fragment emission will be given later. 



B. Light particle emission 

In order to elucidate reaction mechanisms, the emis- 
sion of light particles (Z < 2) is studied. The reactions 
with -"^^^Au are examined first. In Fig. 1201 the emission 
rates of nucleons emitted as a light particle are plotted 
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FIG. 30: Nuclcon emission rate for ^^'^Au is plotted as a funetion of 
time, (a) 47A MeV with NN^m (b) 47A MeV with NNe„p (c) 35A 
MeV with NNiM (d) 26A MeV with NN^a^. Nucleons emitted as 
light particles (Z < 2) are selected. Open circles indicate nucleons 
from the target and solid dots indicate those from the projectile. 



0.2 
0.15 

0.1 
0.05 


0.2 
0.15 

0.1 
0.05 




u 
E 



5 







(a)NNLM 


(b)NNLM 


o ••• 

- • " - 


o^oo 


L*M^ 1 < 




-(c)NN,„p 


- (d)NN,„p 


_ * -.. _ 


"cP 


8 

•M* 1 < 1 


8 



100 200 100 200 300 

Time (fm/c) 



FIG. 31: Similar plots to Fig. | 
(right) at 47A MeV. 



but for 92Mo (left) and ^^Ni 



as a function of time separately for nucleons from the 
projectile and for those from the target. At t ~30 fm/c, 
which corresponds to the full overlap time of the projec- 
tile and the target, light particles start to emerge and the 
emission rate shows a bump at t ~ 80 fm/c for all inci- 
dent energies. At 47A MeV, on the left side of Fig. 1201 the 
emission rate at this bump remains almost equal for nu- 
cleons from the projectile and those from the target. This 
indicates that these nucleons are emitted from a source 
which consists of equal numbers of the projectile nucle- 
ons and the target nucleons, that is, the overlap region 
of the projectile and target. As seen in Fig. [201 (a) and 
(b), the emission rates stay more or less the same when 
the in- medium NN cross section is changed from NN^m 
to NNemp. At the relative NN energy of 50 MeV, the Li- 
Machaliedt formulae gives about 50% larger average NN 
cross section than the empirical formula. This increase 
of the NN cross section enhances the NN collisions but 
also reduces the mean free path of the scattered nucleons. 
The above observation, therefore, suggests that these two 
effects are more or less balanced and the emission rate 
stays more or less same at 47A MeV. As discussed later, 
phenomenologically these nucleons can be described by 
emission from a moving source with half the beam veloc- 
ity, which is experimentally identified as the intermediate 
velocity source. 

At lower energies a similar trend is observed, though 
the emission rates of the target nucleons in the bump 
are about 50% larger at 35A MeV and 20% larger at 
26 A MeV than that of the projectile nucleons. These 
differences are still smaller than the 3:1 ratio of the target 
mass to the projectile mass. The differences may result 
from the fact that the absolute emission rate at these 
energies becomes much smaller than that at 47A MeV 
and it significantly depends on the details of the emission 
process, such as the location of the overlap region in the 
composite system. 

A similar trend is also observed in the other reaction 
systems. In Fig. EH results for ^^Ni and ^^Mo at 47A 
MeV are shown both for NN^m and NNemp- AH emission 
rates show a bump at t ~80 fm/c for nucleons from the 
projectile and those from the target. The strength of 
the emission rate at the bump is ~1.2 both for NN^m 
and NNemp- This value is very similar to that for ^^^Au. 
This is consistent with the above conclusion that these 
preequilibrium particles originate from the overlap zone. 

After the emission of these preequilibrium particles, as 
seen in Fig. OOl the emission rate of nucleons from the 
target increases much faster than that from the projec- 
tile. For the reaction at 47A MeV, the ratio between the 
two reaches a ratio of 3 to 1, essentially equal to the ratio 
of the target mass to the projectile mass, at t ~ 150^200 
fm/c. This time becomes t ^250^300 fm/c for the reac- 
tion at 35A MeV and t > 500 fm/c for the reaction at 
26A MeV. 

The energy spectra of the light particles also provide 
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valuable information on the reaction mechanism. In or- 
der to avoid the Coulomb energy complications for par- 
ticle emissions, neutron energy spectra are studied. In 
Fig. 1321 kinetic energies of neutrons are plotted as a func- 
tion of time. Total and transverse kinetic energies are 
plotted on the left. For all incident energies, both en- 
ergies decrease rapidly as time increases up to t ^ 120 
~150 fm/c and decrease very slowly after that. As seen 
in Fig. at 35A MeV, the composite system under- 
goes multifragmentation at t ~ 200 ~ 300 fm/c, but 
the energy spectra change very smoothly around this 
time range. No effect is observed in the energy spec- 
tra. This indicates that the fragments are already formed 
before this time. On the right, the ratio of the two en- 
ergies are plotted. The ratio is above 2 at t < 100 fm/c 
and rapidly decreases to a value of 3/2 at around 120- 
150fm/c, slightly depending on the incident energy. Then 
the ratio becomes more or less constant after that. The 
value of 3/2 indicates that the neutron emission source 
is fully thermalized. 

In Fig|2ni total and transverse kinetic energies of neu- 
trons are compared for reactions with different targets 
at 47A MeV. Shapes of the spectra are very similar to 
each other for all three targets. The energies for the 
lighter targets are slightly higher, reflecting the semi- 
transparency for these reactions. The similarities in en- 
ergies indicate that the preequilibrium emission mech- 
anism for different reactions is similar, even though the 
semi-transparency is prominent for the reactions with the 
lighter targets. 
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FIG. 32: (Left)Average kinetic energy of neutrons in the center of 
mass system as a function of time for ^^^Au at three incident ener- 
gies. The incident energies are indicated in each paneh Total and 
transverse energy are shown by open squares and circles, respec- 
tively. Energies and times are those when neutrons are identified 
for the first time. (Right) Ratios of the total to transverse energies 
of the left figure. 



C. Fragment emission 

One of the interesting observations in the calculated 
results in Fig. [551 is that the excitation energy of the 
second largest fragment is much lower than that of the 
largest fragment. In order to study this in detail, the 
excitation energies and emission times of fragments are 
examined. Results for ^^"^Au at 47A MeV are plotted as 
a function of fragment mass number in Fig. |^ These 
values are evaluated at the time when each fragment is 
identified for the first time. A fragment, identified at 
a certain time step, often breaks up in the later time 
steps. A fragment is considered as a newly identified 
fragment when the fragment appears with a mass number 
difference of more than 4 mass units between the time 
span of lOfm/c. In that case, the end of the time span 
is taken as the emission time of the new fragment. As 
seen in the upper panel, the distribution shows a peak 
around E*/A 2.5 MeV at A - 15. For the hghter 
fragments their excitation energy rapidly decreases with 
decreasing mass number. For the heavier fragments, on 
the other hand, the excitation energy slowly increases 
with increasing mass number and reaches E*/A ~5'--^6 
MeV at A ~ 60, which is similar to that of the largest 



fragments at t ~ 200 fm/c. The lower panel indicates 
that lighter fragments are identified in a broad range of 
time from 100 fm/c to 500 fm/c, peaking around 250 
fm/c, whereas the heavier fragments of F2 are identified 
mainly around t ~ 200 fm/c. 

This trend is essentially the same for the reactions at 
lower incident energies. In Fig.|2Hl the average excitation 
energy of fragments is compared for ^^''Au at all incident 
energies. For all cases the excitation energy shows a sim- 
ilar trend. At 47A MeV the average excitation energy 
of the lightest fragment starts less than 2 MeV/nucleon 
and increases rapidly to about 3 MeV/nucleon Then 
it increases slowly with increasing mass number. At 
26A MeV, the average excitation energies are about 2 
McV/nuclcon lower than those at 47A MeV. 

Similar observations are also made in the other reac- 
tion systems. In Fig. 1361 the average excitation energies 
of fragments are shown for ^*Ni and ^^Mo at all incident 
energies. The excitation energies of fragments at a given 
incident energy are almost identical and very similar to 
that for Au. This indicates that, in the AMD-V cal- 
culations, cold fragment emission is a common feature of 
intermediate heavy ion reactions. 
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FIG. 33: Average kinetic energy of neutrons as a function of time 
at 47A MeV. Different symbols correspond to different targets, as 
indicated in the figure. Upper set of results shows the total kinetic 
energy Etot and lower set for the transverse energy (Ej^). 
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FIG. 35: Average excitation energy of fragments with Z > 3 is 
plotted as a function of the fragment mass number for ^^^Au at 
three incident energies. The incident energy is indicated in each 
panel. 
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FIG. 34: (a) Excitation energy and (b) emission time of fragments 
with A < 100 are plotted as a function of fragment mass number 
for ^^'^Au at 47A MeV. The largest fragment is excluded in these 
plots. 

It is also interesting to see the excitation energy distri- 
bution for different isotopes with a given atomic charge 
Z. The distribution of produced isotopes significantly de- 
pends on the neutron-proton ratio of the system. In 
Fig.|23isotope distributions are compared between ^^^Au 
and ^^Ni at 47A MeV. As one can see, fragments pro- 
duced in the reaction with ^^^Au are distributed much 
more on the neutron rich side than those with ^*Ni. In 
Fig. I2H1 the average excitation energy of these isotopes 



is plotted for ^^^Au. As seen in the figure, the average 
excitation energy of isotopes for a given Z shows only 
a small variation. This variation is typically within less 
than 1 MeV/nucleon. Very neutron rich isotopes, such 
as ^^O and ^''F, have excitation energies similar to those 
near the /3 stability line. 

The excitation energies studied above are the calcu- 
lated internal energies corrected by the binding energy, 
as given by Eq. E| The binding energy varies, depend- 
ing on the mass and charge of fragments. This variation 
becomes significant for lighter fragments. It is, therefore, 
interesting to see the distribution of the internal energy 
without the correction of the binding energy, in order to 
study the energy partition to the fragments. In Fig. 1391 
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FIG. 36: Similar plots to Fig.EHl but for ^^Ni on left and ^^Mo 
on right sides, respectively. 



the calculated internal energies are plotted as a function 
of mass number for ^^^Au at 47A MeV. The distribu- 
tion shows a quite different trend, comparing to those 
seen in the excitation energies in Figs. 1351 and 1361 espe- 
cially for fragments with A < 15. The internal energies of 
fragments are almost constant and even slightly increase 
for the lighter fragments, whereas their excitation ener- 
gies rapidly decrease. This is because these fragments 
are less and less bound with decreasing mass number. 
For heavier fragments with A > 60, the internal energy 
increases slowly with increasing mass number, which is 
also seen in the excitation energy distribution. This ob- 
servation indicates that each fragment carries away more 
or less the same amount of internal energy per nuclcon, 
which is about -4 MeV/nucleon at 47A MeV. This energy 
becomes about -5 MeV/nucleon at 26A MeV. 

The energy spectra of fragments also provide valuable 
information about reaction mechanisms. In Fig. 1401 the 
average kinetic energies of fragments are shown on the 
left column for ^^^Au. The average kinetic energy in- 
creases almost linearly as mass increases. This linear in- 
crease results partially from the increase of the Coulomb 
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FIG. 37: Isotope distributions of fragments for ^*Ni (upper) and 
-'^^''Au (lower). Lines indicate N=2Z. 



energy for larger fragments and partially from a possi- 
ble expansion energy. The Coulomb energy contribution 
should be similar for the three incident energies, because 
the fragment charge distributions are similar as seen in 
Fig-EHl In order to eliminate the Coulomb energy contri- 
bution from the kinetic energy, the kinetic energy at 26A 
MeV is subtracted from those at 47A MeV and 35A MeV. 
The subtracted results are shown on the right column in 
the figure. At 47A MeV they still show a significant lin- 
ear increase as mass increases, whereas at 35A MeV, the 
linear increase becomes less prominent. As a thermal en- 
ergy contribution to the fragment kinetic energy should 
be constant for different mass fragments, this remaining 
linear increase indicates that these fragments have gone 
through a significant expansion process. The difference of 
the expansion energy between 47A MeV and 26A MeV 
is about 0.5A MeV and O.IA MeV between 35A MeV 
and 26A MeV. The sharp drop of the energy difference 
between these two cases suggests that the expansion en- 
ergy at 26A MeV is very small (< 0.1 A MeV). Therefore 
we can conclude that the approximate expansion energy 
is -0.5A MeV at 47A MeV and ~0.1A MeV at 35A MeV. 

The distributions of fragment kinetic energies are quite 
different for ^^Ni and ^^Mo at 47A MeV, reflecting dif- 
ferent reaction mechanisms. In Fig. ^2 the average total 
kinetic and transverse energies of fragments are compared 
for the different targets at 47A MeV. For ^^^Au, the ra- 
tio of the total to the transverse energy is between 1/2 
3/2, which indicates the existence of a significant radial 
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FIG. 38: Average excitation energy distributions of isotopes for 
'^^'^ Au. The lower plot is the same figure as the upper, but expanded 
for the lighter fragments. 
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FIG. 39: Average internal energy distribution of fragments with 
A > 5 for ^^^Au at 47A MeV. The energy is evaluated when the 
fragments are identified at the first time. 
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FIG. 40: (Left) Average kinetic energy of fragments in the center 
of mass system as a function of fragment mass number for ^^^Au. 
Incident energy is indicated in each panel. Energy is evaluated at 
the time when the fragment is identified at the first time. (Right) 
Subtracted kinetic energy of fragments for Ej. at 47A MeV - E^ at 
26A MeV in the upper panel and for Efc at 35A MeV - Efc at 26A 
MeV in the lower one. 



expansion. For ^^Ni and ^^Mo at 47A MeV, the total 
kinetic energy increases more rapidly as fragment mass 
increases. At A ^--^50, the average kinetic energy for the 
^^Ni target becomes almost double of that for the "'^^^Au 
target. On the other hand the transverse energy does not 
increase as the total kinetic energy increases and stays 
around ~10 MeV. This indicates that the fragment dis- 
tribution is stretched along the beam axis. This results 
from the semi-transparency, as reported in the previous 
work 0. 



D. Thermal equilibration time 

In Fig. 1201 for ^^^Au, the ratio of the emission rate 
between the target nucleons and the projectile nucleons 
reaches a factor of 3 at t --200 fm/c at 47A MeV, t --300 
fm/c at 35 A MeV and t > 500 fm/c at 26 A MeV. This 
factor of 3 is roughly the ratio of the target mass to the 
projectile mass and expected for nucleon emissions from a 
thermalized system. However this observation does not 
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FIG. 41: Average total and transverse fragment energies as a func- 
tion of fragment mass number in the center of mass system for 
reactions at 47A MeV with the different targets. The target is 
indicated in each figure. Dots indicate the total energy and open 
circles indicate the transverse energy. 



mean that the system reaches thermal equihbrium un- 
til these late stages, because a significant amount of the 
target nucleons are carried away by cold light fragments 
before these times. These fragments decay very slowly 
and do not contribute to the emission rate at the early 
stages. 

In order to elucidate the thermal equilibration time, 
the ratio of the excitation energies of the largest fragment 
(Pmax) and the second largest fragment (F2) is plotted 
in Fig. 021 When all fragments with Z > 3 are taken into 
account for F2 (circles), the ratio gradually increases and 
reaches a plateau at t --200 fm/c at 47A MeV, t ~250 
fm/c at 35A MeV and t > 500 fm/c at 26A MeV. These 
times are very similar to those for which the nucleon emis- 
sion rates become equal to the 3:1 ratio. On the other 
hand when only fragments with A > 30 are taken into 
account for F2 , the ratio behaves in a drastically different 
way and stays more or less constant from the earliest time 
when these heavier fragments are identified. This time is 
t = 130 fm/c at 47A MeV, 160 fm/c at 35A MeV and 180 
fm/c at 26 A MeV. The values of the ratios are slightly 
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FIG. 42: Ratio of the excitation energy of the second largest frag- 
ment ( F2 ) and that of the largest fragment ( Fmai ) as a function 
of time for ^^^Au. Open circles represent results when all fragments 
with Z > 3 are taken into account as F2 and solid dots indicate 
results when only fragments with A > 30 are taken into account as 
F2. Incident energy is indicated in each figure. 



lower than 1.0, because the excitation energies of A ~ 
30 ~ 40 are still slightly lower than those of the heavier 
fragments, as is seen in Fig. 1351 This observation indi- 
cates that the system is already thermally equilibrated 
at the earliest times that these heavier fragments can be 
identified. 

From the nucleon emission rates in Figs. 1201 and 1311 

one can see that the preequilibrium emission ceases at t 
~120 fm/c at 47A MeV and slightly later at lower inci- 
dent energies. The ratios of the neutron total to trans- 
verse energies on the right side of Fig. |221 reach a value of 
3/2 at these times and remain constant after that. The 
value 3/2 suggests that the neutron emission source is 
fully thermalized. These times are consistent with those 
observed from the earliest emission times of F2 with A > 
30. Therefore from these observations we conclude that 
the excited nucleus, remaining after the preequilibrium 
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emission, reaches thermal equihbrium at t '--'120 fm/c at 
47A MeV, t ~140 at 35A MeV and t -160 fm/c at 26A 
MeV. 

E. Multifragmentation and cold fragment emission 

In this subsection we briefly summarize the observa- 
tions for the calculated central collision events and draw 
a scenario for multifragmentation in Fermi energy heavy 
ion reactions. 

For light particle emission we summarize as follows: 

• The emission rate of nucleons emitted as light par- 
ticles shows a maximum at t ~ 80 fm/c for all re- 
actions studied here. At a given incident energy 
all reactions show similar rates for these preequi- 
librium particles emitted from the projectile and 
from the target. The emission rates are insensitive 
to the change of the in-medium NN cross section. 
The absolute emission rate decreases significantly 
as the incident energy decreases, but even at 26A 
MeV a similar rate is observed for the projectile 
and target nucleons. These facts suggest that the 
preequilibrium light particles are emitted from the 
overlap zone of the projectile and the target. 

For fragment emission we summarize as follows: 

• Charge distributions of fragments with Z < 20 are 
very similar to each other in the experiments and 
in the calculations, regardless of the target. 

• Cold fragment emission is generally observed for all 
reactions studied here. The excitation energy for A 
< 10 is E*/A = 1—2 MeV. It increases linearly as 
mass increases and reaches to that of the largest 
fragment at A ~ 60, remaining essentially constant 
for heavier fragments. The rapid increase of the 
excitation energy for fragments with A < 15 origi- 
nates from the variation of the binding energy. The 
internal energies of fragments show a flat distribu- 
tion for A < 60. 

• For ^^^Au a significant radial expansion is ob- 
served at 47A MeV. The expansion energy de- 
creases quickly as the incident energy decreases. 
The approximate expansion energy is ~0.5A MeV 
at 47A MeV and -0.1 A MeV at 35 A MeV. 

• For the reactions with '^^Ni and ^^Mo at 47A MeV, 
a significant semi-transparency is observed. 

For the thermal equilibration time we summarize as 
follows: 

• The thermal equilibrium of the system is estab- 
lished around t ~ 120 fm/c at 47A MeV, t - 140 
fm/c at 35A MeV and t ~ 160 fm/c at 26A MeV. 



From these observations we can draw the following sce- 
nario for the multifragmentation in the Fermi energy do- 
main: 

(1) The projectile enters into the target nucleus and 
creates a hot overlap zone. 

(2) The overlap zone decays quickly by emitting 
fast nucleons before thermal equilibrium with the 
surrounding target nucleons is established. 

(3) During the preequilibrium emission, the re- 
maining system starts to equilibrate. Nucleons 
close to each other in the phase space start to form 
cold fragments which coexist with a hotter nucleon 
gas. In this stage, cold fragments share almost an 
equal internal energy per nucleon with each other. 

(4) The system continues to expand and undergoes 
multifragmentation with cold fragment emission. 
This multifragmentation process is accelerated by 
additional dynamical processes, such as the semi- 
transparency for ^^Ni and ^^Mo at 47 A MeV and 
the radial expansion for ^^^Au at 47A MeV. These 
additional dynamical processes drastically change 
the kinetic properties of the emitted fragments. 

VII. DISCUSSION 

For the above multifragmentation scenario, the exis- 
tence of a hot overlap zone and cold fragment formation 
play important roles. For the existence of the hot overlap 
zone, many experimental evidences have been accumu- 
lated. Since the pioneering work of Awes et at, moving 
source analyses have been applied to many reaction sys- 
tems at intermediate energies |43,|43,|5fl|5lllS2; 53]. The 
existence of the intermediate velocity source has been 
commonly observed in these analyses. The source veloc- 
ity extracted from these analyses is about half the beam 
velocity, regardless of the target mass. This has been in- 
terpreted as indicating that the source consists of equal 
amounts of nucleons from the projectile and the target. 
The apparent temperature of the source increases linearly 
as the incident energy increases, and is more or less in- 
dependent of the target mass and the centrality of the 
reaction These observations are consistent 

with preequilibrium emission of light particles from the 
overlap zone, discussed in the previous section. 

A more direct indication for the preequilibrium parti- 
cle emissions from the overlap zone was reported recently 
by Verde et al. js^l . In that work, two-proton correlation 
functions were studied using a newly developed imaging 
source technique. The source image of two protons was 
derived numerically from the two proton correlation func- 
tion, without assuming a source shape, such as a Gaus- 
sian distribution. This method is capable of determining 
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the source size in the environment in which fast (pree- 
quilibrium) and slow (evaporated) components coexist in 
spectra. The model was applied to ^"'N + ^^^Au at 35A 
MeV and the source size of the preequilibrium protons 
was extracted. The extracted source size is Ri ~ 3 fm 

2 

for all proton energy ranges. (Ri is the radius at a half 
density). This source size is comparable to the size of the 
projectile and much smaller than the size of the target. 

Cold IMF emission has been reported in many heavy 
ion collisions. The temperature of IMFs has been deter- 
mined from the population of the excited states [s^ 
I'^ these studies temperatures of T ~3~5.5 MeV 
have been typically obtained in a variety of intermedi- 
ate heavy ion reactions. If one uses the relation E*=aT^ 
with a=A/10, these temperatures lead to fragment ex- 
citation energy of E*/A = 1~3 MeV. Marie et al. also 
reported cold fragment emissions observed in LCP (light 
charged particle)-IMF correlation studies |H3 . The aver- 
age charged particle multiplicities for each emitted IMF 
were evaluated in the ^^^Xe + Sn reaction at 50A MeV. 
Using a statistical cascade code, the average excitation 
energy was determined for each IMF. When the neutron 
to proton ratio of IMFs is assumed to that of the system, 
the excitation energy of E*/A ~ 3 MeV is obtained for 
fragments with ?> < Z < 20. This energy is significantly 
lower than the excitation energy ( E*/A ~ 12 MeV) of 
the system, but consistent with the value calculated in 
this study. The correlation experiment has recently been 
extended from 25A MeV to 150A MeV in the same re- 
action system The extracted excitation energy of 
fragments is E*/A ~ 3 MeV, independent of the incident 
energy. 

In each of the studies in Refs. [s^ and [H^l, however, 
the extracted excitation energies are slightly decreasing 
as Z increases, whereas the results of the calculations in 
Figs.l35landl36lshow a clear increase of the excitation en- 
ergy as fragment mass increases, especially for the lighter 
fragments. This difference may result from the fact that, 
in both of the experimental studies, the fragment isotope 
distribution was neglected and neutron emission from the 
fragment was not measured. Thus the average primary 
fragment mass for a given average primary fragment Z 
was calculated by (l-l-f)Z, where f is N/Z of the frag- 
ment and it was assumed in Refs. [HI and [HOl that N/Z 
is that of the composite system. If N/Z is assumed to 
be that of the valley of the stability rather than that 
of the composite system, the extracted excitation ener- 
gies become less than 1.5 MeV/nucleon for most of the 
fragments. As seen in Fig. |^ for a given Z the calcu- 
lated distributions of fragment mass are broad and the 
N/Z value at the peak of the mass distribution is slightly 
less than that of the system. For example, the distribu- 
tion of Oxygen isotopes has a peak at ^*0, or N/Z = 
1.25, for ^"^Zn -I- ^^^Au, which has a composite system of 
N/Z = 1.39. Since the experimental values are sensitive 



to the determination of the primary fragment mass and 
hence the number of emitted neutrons, one needs to take 
into account the isotope distribution of the fragments for 
the evaluation of the average excitation energy in order 
to make more detailed comparisons between the experi- 
ments and the calculations. 

Recently Cussol suggested that the limitation of the 
excitation energy of fragments may be related to the neu- 
tron (or proton) separation energies of the fragments, us- 
ing a Classical Molecular Dynamics model |61|. In most 
heavy ion reactions, the composite system is neutron rich. 
When light fragments are formed in an equilibrated sys- 
tem, the neutron/proton ratio of light fragments tends to 
be far away from the (3 stability line and these fragments 
have small neutron separation energies. These fragments 
can not hold excitation energy higher than the neutron 
separation energy and therefore they limit the average ex- 
citation energy of the fragments. In AMD-V simulations, 
however, no systematic correlation is observed between 
the neutron/proton ratio and the excitation energy, as 
shown in Fig.|2Hl The calculated results also show no sig- 
nificant difference between the neutron rich system (^^Zn 
+ ^^^Au ) and more or less symmetric system (^''Zn -|- 
58Ni). 

The coexistence of cold fragments and gaseous nucle- 
ons has recently been suggested by Campi et al., using 
a Classical Molecular Dynamics (CMD) model In 
their CMD model, particles interact through a Lennard- 
Jones potential, which is given by 

V{r)=Ae{{-f^-{-f}, (19) 
r r 

where e and a are energy and length scales, respectively. 
In the reference a thermalized system is provided with a 
given number of particles and a fixed excitation energy 
within a small container in high density. At time=0 the 
container is removed and the system is allowed to expand 
freely. It undergoes multifragmentation. Fragment size 
and temperature are examined during the process. Since 
the particles interact through the above interaction, they 
tend to move at a distance ro from each other ( V(r) be- 
comes minimum at r=ro ) . Using proximity in the phase 
space and separating gaseous nucleons by the Hill cri- 
terium (potential energy -I- kinetic energy > 0.), clusters 
are identified in the gaseous particle environment. By 
examining the temperature of gaseous particles and the 
internal temperatures of the clusters during the expan- 
sion, it is revealed that the temperatures are much lower 
than the temperature of the gaseous particles. The dis- 
tribution of cluster size is almost independent of the time 
of the expansion, indicating that the cluster size is inde- 
pendently determined of the particle density. 

In AMD-V, nucleons are moving freely in a mean field, 
but anti-symmetrization tends to maintain clusters with 
nucleons which are close each other in the phase space. 
This is illustrated in Fig. 021 In the figure the root-mean- 
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as radial expansion and semi-transparency, drastically 
change the kinematic characteristics of the emitted frag- 
ments. A drastic change in kinematic distributions has 
been reported in ^"^Ar induced reactions [g^I- In that ex- 
periment '*°Ar beam was bombarded on Cu, Ag and Au 
targets over a wide range of incident energy (8 A -115A 
MeV) . For the central collisions below 44A MeV, the ob- 
served nature of the reactions indicates a fusion-like pro- 
cess, whereas a multi-body spray of IMFs is observed at 
higher incident energies. This observation is consistent 
with the present results for ^*Ni. At 26A MeV the pro- 
jectile fully stops in the target and fragment emissions 
become more or less isotropic. On the other hand at 47A 
MeV, a significant semi-transparency spreads the frag- 
ments at forward and backward directions. 

The existence of a significant expansion has been 
suggested based upon the determination of the aver- 
age kinetic energy of IMFs by the INDRA collabora- 
tion [H^ . Issf . The extracted energy of the possible ex- 
pansion is about 0.5A MeV for Xe-I-Sn at 50A MeV, con- 
sistent with the observation in Fig.QHlat 47A MeV. How- 
ever the value in the above references has been evaluated 
based on the statistical multifragmentation model, which 
is not supported in the present work. 



VIII. SUMMARY 



FIG. 43: Time evolution of the average root-mean-square radius 
of 14 nucleons, which end up ^^C at t=480 fm/c for ^^^Au. See 
details in the text. 

square radius ( Krms ) of "'^^C, which is the most abundant 
light fragment for ^^^Au, is examined. 14 nucleons in 
""^^C are identified at t=480 fm/c. For these identified 14 
nucleons, Rrms is calculated as a function of time. At 
t < 0, the time before the two nuclei touch each other, 
Rrms shows large values because about 3^4 nucleons on 
average originate from the projectile and the rest from 
the target. At 47A MeV, Krms reaches 2.6 fm at t = 30 
fm/c, the time of the full overlap, and after that, R^ms 
decreases very slowly. The change in R^ms between t 
= 30 fm/c and t = 480 fm/c is ~10%. This trend is 
essentially the same for lower incident energies. At 26A 
MeV, Krms becomes 2.4 fm at t = 50 fm/c and 2.2 fm 
at 480 fm/c. The slightly smaller Rrms, compared to 
that at 47A MeV, reflects the fragment lower excitation 
energy at 26 A MeV. This indicates that nucleons in light 
fragments are already close together at the time of the 
overlap and move together until they are identified as 
an isolated fragment. The excitation energies of light 
fragments are small, independently of the emission time, 
as seen in Fig. I34r b'l. This indicates that the fragments 
stay cold and coexist with the hot nucleon gas. 

Different dynamics involved in the reactions, such 



A detailed study has been presented for three reaction 
systems, ^'^Zn + ^^m, ^^Mo and ^^^Au, at three incident 
energies, 26A, 35A, 47A MeV. Multiphcity distributions, 
charge distributions, energy spectra and velocity distri- 
butions of the reaction products have been measured. 
Detailed comparisons have been made between the exper- 
imental results and those of the calculations with two dif- 
ferent effective interactions, corresponding to a soft EOS 
(K=228 MeV) and a stiff EOS (K=360 MeV). For light 
particles the main characteristic features of the above 
observables are generally well described by all three cal- 
culations, although some significant discrepancies are ob- 
served in all cases. Experimental proton multiplicity dis- 
tributions, for example, favor the calculations with the 
soft EOS, whereas those of the a particles favor the cal- 
culations with the stiff EOS. For the fragments with Z > 
3 the experimental energy spectra clearly favor the cal- 
culations with the stiff EOS, although the multiplicity of 
IMF's is overpredicted by a factor of 1.5-2 for ^^Ni and 
^^Mo and 2-3 for ^^^Au. The velocity distributions of a 
particles and IMF's also show a clear difference between 
calculations with the soft EOS and the stiff EOS, relating 
to the velocity distribution of nucleons in an early stage 
of the reaction. The experimental velocity distributions 
of these particles clearly favor the calculations with the 
stiff EOS. Two different formulations for the in-medium 
NN cross sections were used in this study, but the effect 
of changing the NN cross section is rather small on the 
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observables and no conclusive evidence favoring one of 
these formulations has been found. 

The mechanism of the multifragmentation process was 
explored, using central collision events of AMD-V, cal- 
culated with the stiff EOS. Preequilibrium light particle 
emission from the overlap zone and multifragmentation 
with cold fragment emission are commonly observed in 
all reactions. The thermal equilibration times observed 
in the calculated results are t ~ 120 fm/c at 47A MeV, t 
~ 140 fm/c at 35A MeV and t - 160 fm/c at 26A MeV. 
Cold fragments are formed by nucleons close to each 
other in phase space at early stages and they stay cold 
in a hotter nucleon gas, exhibiting an almost equal inter- 
nal energy per nucleon. Then the system expands and 
undergoes multifragmentation with cold fragment emis- 
sion. For ^^"^Au at 47A MeV a significant radial expan- 
sion takes place and for ^*Ni and ^^Mo at 47A MeV, semi- 
transparency becomes significant. The kinematic charac- 
teristics of emitted fragments change drastically, depend- 
ing on the additional dynamics involved. Many existing 
experimental results are consistent with the mechanisms 
suggested by the calculated AMD-V events. 
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